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ABSTRACT 

We carried out a comprehensive far- ultraviolet (UV) survey of ^^CO and H2 column densities along 
diffuse molecular Galactic sight lines in order to explore in detail the relationship between CO and H2. 
For this survey we measured new CO abundances from absorption bands detected in Hubble Space 
Telescope spectra for 62 sight lines, and new H2 abundances from absorption bands in Far Ultraviolet 
Spectroscopy Explorer data for 58 sight lines. In addition, high-resolution optical data were obtained 
at the McDonald and European Southern Observatories, yielding new abundances for CH, CH+, and 
CN along 42 sight lines to aid in interpreting the CO results. A plot of log 7V(C0) versus log iV(H2) 
shows that two power-law relationships are needed for a good fit of the entire sample, with a break 
located at log N{GO, cm^^) = 14.1 and log iV(H2) — 20.4, corresponding to a change in production 
route for CO in higher-density gas. Similar logarithmic plots among all five diatomic molecules allow 
us to probe their relationships, revealing additional examples of dual slopes in the cases of CO versus 
CH (break at log N = 14.1, 13.0), CH+ versus H2 (13.1, 20.3), and CH+ versus CO (13.2, 14.1). 
These breaks are all in excellent agreement with each other, confirming the break in the CO versus 
H2 relationship, as well as the one-to-one correspondence between CH and H2 abundances. Our new 
sight lines were selected according to detectable amounts of CO in their spectra and they provide 
information on both lower-density (< 100 cm~'^) and higher-density diffuse clouds. The CO versus 
H2 correlation and its intrinsic width are shown to be empirically related to the changing total gas 
density among the sight lines of the sample. We employ both analytical and numerical chemical 
schemes in order to derive details of the molecular environments. In the denser gas, where C2 and CN 
molecules also reside, reactions involving C"'' and OH are the dominant factor leading to CO formation 
via equilibrium chemistry. In the low-density gas, where equilibrium-chemistry studies have failed to 
reproduce the abundance of CH"*" , our numerical analysis shows that nonequilibrium chemistry must 
be employed for correctly predicting the abundances of both CH+ and CO. 
Subject headings: astrochemistry — ISM: abundances — ISM: molecules — ultraviolet: ISM 



1. INTRODUCTION 

The most abundant molecule in the Cosmos, H2, has 
no permanent dipole moment and is thus lacking permit- 
ted pure rotation and vibration-rotation transitions. On 
the other hand, rotational transitions of CO, such as J = 
1—0 at 115 GHz, have been routinely and extensively ob- 
served in molecular clouds, whi ch are also too cold for de- 
tection of excited levels of H2 (jPapadopoulos et"al]l2002l 
and references therein). Radio telescopes are thus used 
to map CO emission in the interstellar medium (ISM) in 
order to delineate global distributions of molecular clouds 
in our Galaxy and in other galaxies. 

It is an empirical and theoretical foundation of radio 
mapping that the velocity-integrated emission intensity 
of CO (Wco, in units of K km s~^) from a molecular 
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cloud is proportional to the total virialize d mass of the 
cloud, and hence to its hydrogen content (jLarsonl[l981l : 
lYoung fc Scovilldll991h . Radio astronomers thus utilize 
Wco as a proxy for H2, by employing the "X-factor" 
Xco = A^(H2)/VFco(l - 0), where N is the observed 
column density and Xqq will be assumed hereafter to be 
in units of 10^*^ cm^^ (K km s~-^)~^. 

An average value is ascribed to gian t molecular clouds 
(GM Cs) in the Milky Way , X^p » 4 ("Young & ScOTmi 
1982: JDickman et al.lflQSl . iPolk et a l. (1988) found sig- 
nificant millimeter-wave emission from CO that was not 
associated with GMCs. This lowers the value of Xqq, 
having a mean value in the solar neighborhood of 1.8 
± 0.3 (jDame et al.ll200ll ). A dependence of ATco on the 
metallicity (primarily C/H) of the gas has been found in, 
e.g., studies of 7-ray emission across the Galaxy. This 
intensity is the product of interactions between cosmic 
rays and "stationary" gas; thus the 7-ray intensity is 
proportional t o the amount of gas along the line of sight. 
IStrong et al.l ()2004l ) found that Xco varies between 0.4 
and 10.0 from the inner Galaxy to its outer regions, in- 
dicating lower gas metallicity in the outer Galaxy. 

The GMCs, as well as the smaller dark clouds, are 
opaque enough that the CO in their cores is not disso- 
ciated by the interstellar far-ultraviolet (far-UV) radia- 
tion field. Consequently, CO is unobservable in the UV, 
necessitating its detection via millimeter-wave emission. 
On the other hand, diffuse molecular clouds (as well as 
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envelopes of dark clouds) have visual extinctions lower 
than 5 mag, enabling a direct determination of column 
densities via line absorption in the UV. Both CO and 
H2 are photodissociated by far-UV radiation, resulting 
in a variable Xqo that depends on the efficiency of self 
shielding (as well as mutual s hielding) of the two species 
(|van Dishoeck fc BlaclJll988[ ). and thus on their column 
densities. One of our goals here is to study the behavior 
of Xc n under diffuse-I S M co nditions. 

The iFeder man et al.l (1980) study of diffuse interstel- 
lar clouds showed that there is an approximate quadratic 
relationship between CO and H2, such that A^(CO) oc 
[JVfHg )]^, with B « 2. Furthermore, IFederman et all 
(|198Clf ) remarked that for a group of lower- A'' sight lines, a 
shallower slope with B w 1.5 was more appropriate, thus 
signaling the possibility of B varying with N. Indeed, for 
sight lines with log 7V(C0) ^ 16, or values higher than 
those that were available to IFederman et al.l 1)19801 ) , an 
even ste eper relati o nship with B kj 3 was found in the 
studv of|Pa"n et all (l2005f)^ Both[BlI7gh et all |2007) and 
ISonnentrucker et al.l ( 20071 ) confirmed the steepening of 
the slope near log iV(CO) « 15, attributing it to self- 
shielding of CO. Thus a variable power law in CO ver- 
sus H2 shows that the abundance of CO relative to H2, 
and hence Xco, are increasing with A^(H2) along dif- 
fuse molecular sight lines. In this paper we explore the 
CO versus H2 relationship in more detail, namely, the 
trend of CO (as well as of CH, CH+, and CN) versus H2 
in diffuse molecular clouds and the dependence of these 
correlations on A'^(H2) and on physical parameters, such 
as the total gas (hydrogen) density, nn- 

Our study differs from the recen t work of iBurgh et all 
(j2007( ) and ISonnentrucker et all (I2007D in a num- 
ber of ways. Fi rst, unlike iBurgh et aLJ and like 
ISonnentrucker et al.l . our methodology is based on pro- 
file fitting of the data, with detailed decomposi tion into 
cloud component structures (jSheffer et al.ll2007l ). We do 
not consider apparent optical depth nor curve of growth 
treatments, which are alway s less preferred to spectrum 
synthesis by profile fitting (So nnentrucker et al.ll2007l ). 
nor do we build a grid of mod els to look for a so lution 
with a single effective 6- value (jBurgh et al.ll2007h . Sec- 
ond, we follow up on these measurements of N values 
with two methods of chemical analysis, analytical and 
numerical, in order to derive riu at the sites where CO 
is detected. Third, we are able to discern two regimes 
of CO formation in terms of rifj. CO is associated with 
the similarly heavy diatomic mole cules C2 and CN inside 
denser and colder clumps of gas |Federman et al.lll994l : 
iPan et aP 120051 : ISonnentrucker et "all l2007f ) , whereas in 
low-density clouds CO is related to the form ation and 
chemistry of CH+ (|Zsarg6 fc FedermanI 120031 ). Remark- 
ably, this transition in the photochemistry of CO will be 
shown (§ 3) to affect also the trends of other correlations 
among the diatomic molecules analyzed here. 

In § 2 we detail our sources for data and our methods of 
reduction and analysis. Next, in § 3 observational results 
are presented in terms of derived component structures 
and correlations between molecular column densities. In 
§ 4, 5, and 6 we explore some of the physical conditions 
of the CO-harboring gas in terms of empirical relation- 
ships, analytical chemical analysis, and detailed numeri- 
cal modeling with Cloudy, respectively. A discussion will 



be given in § 7, followed by the conclusions in § 8. 
2. DATA AND MODELING 

Our primary effort was to detect and measure A^(CO) 
for new sight lines from archival Hubble Space Telescope 
(HST) data. For most of these sight lines the value of 
A^(Il2) was already kn own from previous s urveys with the 
Copernicus satellite () Savage et al.l Il977l ) or Far Ultra- 
violet Spectroscopic Explorer (FUSE) (JRachfor d et al. 
12002"; 'Andr e et al.ll2003l : ICartledge efadmO'k iPan et al. 
.2005). However, for consistency, we determined A^(H2) 
also for sight lines with previously published results. 
Only two new sight fines (HD 36841 and HD 43818) lack 
any A^(H2) data. For these, predicted values of A^(Il2) 
will be provided in § 4.4 after exploring the II2 relation- 
ships with CO and CH. We obtained new high-resolution 
optical spectra of CH+, CH, and CN by observing 42 
sight lines at either McDonald Observatory or at the Eu- 
ropean Southern Observatory (ESO). Results for ^'^CO, 
which is also present in the HST spectra, were published 
in ShcflPer et al. (2007). 

Table 1 provides a list of all sight lines in terms of stars 
observed, their spectral types, visual magnitudes. Galac- 
tic coordinates, local standard of rest (LSR) corrections, 
E{B — V) reddening values, and heliocentric distances. 
Table 2 lists the UV data sets from HST and FUSE for 
our stellar targets, and STIS optical setups in terms of 
gratings and apertures. 

2.1. HST Data 

Initially, our sample included 66 sight lines with- 
out previous measurements of A^(CO). The UV data 
for 63 of these consist of archival STIS observations, 
from which we extracted spectra oi A — X bands of 
CO between 1229-1544 A. The remaining three sight 
lines have archival GHRS data. Results on A^(^^CO) 
were subsequently publishe d for 23 sight lines: 12 in 
IBurgh et all (|2007D . three in ISonnentrucker et al.l (I2007D 
(with one sight line in com mon with IBurgh et aLh , and 
12 in 'Sheffer et al.l (|2007l ) (with two in c ommon with 
[Burgh et al. a nd one in c ommon with both iBurgh et al.l 
and Sonne ntrucker et al.l ). Thus, this paper presents new 
N(^'^CO) results for 43 sight lines. To the entire CO 
sample of 66 sight lines we added previously published 
A^(CO) values for 48 directions, yielding a sample of 114 
sight lines with UV data. Figure 1 presents the view 
of CO absorption along two HST /STIS sight lines that 
differ by a factor of 700 in N{CO). 

2.2. FUSE Data 

Our initial sample of 58 sight lines was obtained from 
archival FUSE observations of H2 absorption at A < 1100 
A. Of these, 33 sight fines did not have published A^(H2) 
results. In the meant ime, JV(H2) results were published 
for five sight lines in 'Burgh et al.l (|2007t ) and five more 
in Sheffer et al. (2007) (with a single sight line in com- 
mon with Bur gh et al.l ). This paper, therefore, presents 
first A^(H2) results for 24 sight lines. As described in 
IFederman et al.l ()2005f ). our iV(H2) values are obtained 
from spectrum synthesis of the (2—0), (3—0), and (4—0) 
bands of the Lyman B — X transitions of H2. The to- 
tal column density A^(H2) listed is based on the absorp- 
tion from all rotational levels with J" = through 5. 
Roughly 95% of the total is found in the two J" ~ and 
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1 ground states of para- and ortho-H2, respectively. Fig- 
ure 2 presents a sample of two FUSE sight lines with H2 
absorption profiles that differ by a factor of 13 in iV(H2). 
The spectral coverage of F USE also c o ntains 
absorption features from CO (jSheffer et al.1 120031 : 
iCrennv &: FedermanI l2004f ) . Thus for three sight lines 
(HD 208905, HD 209481, and HD 209975) with no HST 
spectroscopy we determined A^(CO) from the B ~ X 
(0-0), C - X (0-0), and E - X (0-0) bands of CO, 
as well as confirmed the CO content along the line of 
sight towar d HD 200775 that was previously based on 
lUE data (jKnauth et al.ll2001h . The former three stars 
were included in t he high-resolution optic al study of CH, 
CH+, and CN bv iPan et alJ ((200I [20051) . 

2.3. McDonald Data 

High- resolution optical observations of CH+, CH, and 
CN were obtained with the 2dcoud e cross-dispersed 
echelle spectrometer ()Tull et al.l I1995D for the purpose 
of deriving cloud structure templates for CO and H2, 
without which one cannot reliably derive the line opti- 
cal depth for sight lines with very high column densi- 
ties. In addition, the echelle spectra included absorption 
from Ca II and CN, the first provides a high signal-to- 
noise confirmation of the cloud structure, while 7V(CN) 
is used to model the total gas density in the absorbing 
cloud, based on the CH and C2 chemical reaction network 
described in § 5.1. Additional absorption from CH+ pro- 
vides a check on the component structure for directions 
with low molecular concentrations. 

Sight lines toward 20 stars were observed at i? ~ 
170,000 with the 2.7-m Harlan J. Smith telescope at 
McDonald Observatory, Texas, during observing runs in 
2004 January, October, and December, and 2005 May 
and October. Each echelle exposure included nine or- 
ders, which simultaneously recorded two atomic transi- 
tions, Ca I at 4226 A and the K hue of Ca II at 3933 A, as 
well as absorption lines from three molecules, CH at 4300 
A, CH+ at 4232 A, and CN at 3784 A. Two-dimensional 
reduction tasks in IRAF were used to correct these for 
bias, scattered light, pixel-to-pixel variations, and finally 
to calibrate the wavelength scale based on accompany- 
ing exposures of a Th-Ar lamp. The latter step yielded 
residuals smaller than 0.001 A, or < 0.07 km s~^. 

2.4. ESO Data 

For the 17 stars in our sample that are located too 
far south to be observable from McDonald, data were 
obtained at ESO^ in Chile. Five more sight lines were 
added to the ES Q observing prog r am to complement CH 
results given in 'A ndersson et al.l ()2002[ ) with new CH+ 
and CN acquisitions. For the 22 sight lines we obtained 
exposures on CH for 16 sight lines, on CH+ for 19 sight 
lines, and on CN for five sight lines. The observations 
were carried out at the 3.6-m telescope at LaSilla in 2005 
June and 200 6 June using the Coude Echelle Spectro- 
graph (CES) (lEnardlll982(l . The CES is fed with an op- 
tical fiber with an aperture of 2" on the sky and provides 
R of 220,000. The data were bias subtracted, flat-fielded, 
and rebinned to a linear wavelength scale using the MI- 

^ Based on observations collected at the European Southern Ob- 
servatory, La Silla, Chile, programs nos. 075.C-0025(A) and 077. C- 
0116(A). 



DAS long package. Figures 3, 4, and 5, provide com- 
parisons of sight lines with weak and strong absorption 
from CH+, CH, and CN, respectively. HD 23478 and HD 
210121 were observed at McDonald, whereas the acqui- 
sition of data for HD 99872 and HD 116852 occurred at 
ESO. 

2.5. Ismod.f Spectrum Synthesis 

We used the Y.S. code, Ismod.f, to model Voigt ab- 
sorption profiles via spectrum synthesis and automatic 
rms-minimizations of (data minus fit) residuals. Be- 
sides presenting the data. Figs. 1-5 all include spectrum 
synthesis fits that were performed with Ismod.f. The 
basic absorption equa t ions w ere adapted in 1990 from 
iBlack fc van Dishoeckl (|1988[ ). Besides fitting radial ve- 
locity, excitation temperature (Tox) , and total N for any 
absorption feature, Ismod.f provided solutions for cloud 
structures along each sight line, i.e., the number of cloud 
components, their relative shifts and fractions, and their 
Doppler widths (&- values). This information is critical for 
proper evaluation of large optical depth effects and has 
to be derived ab initio whenever not known from previ- 
ous investigations. Table 3 presents all cloud components 
that we were able to identify via molecular absorption, 
while Figs. 3 and 5 present spectra of CH+ and CN to- 
ward HD 23478, with each species clearly showing two 
cloud components. Our criteria for detection were a 2-a 
limit for molecular col umn de nsity, as well as simultane- 
ous detection in Ca II (jPan et al. 2005) . All radial veloc- 
ities have been transformed from the heliocentric scale 
to the LSR reference fra me. For CO transition st r ength s 
we used the /-values of IChan. Cooper, fc BrionI (|1993f l. 
which in a global sense have been verified to a level of 
a few percent by Eidelsbcrg ct al. ( 1999 ,), whil e wave - 
lengths were obtained from iMorton fc Noreaul (|1994( l. 
Both /-values and wa velengths for H2 were obtained from 
lAbgrall et all (|1993allbl ) via files available on the Dr. Mc- 
Candliss website*. As for the species with optical transi- 
tions, the corresponding input values for Ca I, Ca II, CH , 
CH+, and CN were taken from Table 3.4 in Pa3 (2002). 

A subset sample of i V(^^CO) and A^("CO ) for 25 sight 
lines was published bv lSheffer et al.l (|2007| ). who showed 
that a minimal number of absorption bands is needed for 
a robust modeling of A^(CO). Specifically, a few bands 
that are optically thin and a few that are optically thick 
should be simultaneously synthesized to yield a good 
measure of N and of those parameters that affect line 
saturation in the bands, such as b and Tox- Based on 
that sample of CO results, we find that our uncertain- 
ties in N range from a few to ~20%. Thus to be on the 
conservative side, we shall assume that the 1-ct errors are 
±20%, plotting this value in all of our relevant figures. 

All H2 lines from the J" < 5 levels were modeled 
simultaneously with Ismod.f. The major difference in 
modeling methodology between CO and H2 is that for 
the latter we do not attempt to fit any parameters that 
are associated with the cloud component structure along 
the line of sight. This is a direct result of the relatively 
low spectral resolution of FUSE {R - 20,000) as well 
as the restricted available range (compared to CO) in /- 
values for the Lyma n and Werner ba nds of H2. Since it 
had been discovered (jFedermanll 1 9821 ) and will be verified 

8 www. pha.jhu.cdu/'~stcphan/h2ools2. html 
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in § 3.3 that 7V(CH) and iV(H2) have a Unear relation- 
ship, our method is to apply any cloud structure already 
known from high-resolution CH data to the modeling of 
H2, while keeping such structure para meters fixed dur- 
ing the fit. As in ou r previous work (Pan et all l2005l : 
iFederman et al.ll2005( l we prefer known cloud structures 
to effective 6- values as the proper solution to H2 line sat- 
uration. This is in contrast to CO, where parameters 
such as relative strength and width of components are 
allowed to vary during th e fit. Further det ails can be 
found in our earlier paper (jSheffer et al.ll2007[ ) . In a sim- 
ilar fashion to CO, we shall show global ±20% l-cr error 
bars in all our figures that present values for 7V(H2). This 
uncertainty is consistent with published results. 

Whenever repeated multiple exposures are available, 
we combined them in wavelength space for an improved 
signal-to-noise (S/N) ratio. In addition, when a fea- 
ture (band) appears in two adjacent orders, we combined 
them after correcting for any small wavelength inconsis- 
tencies by subtracting the wavelength shifts as measured 
from the absorption line positions. All our reductions 
were in IRAF and STSDAS. A single Gaussian was used 
to describe the instrumental profile of STIS, but R was 
allowed to be a free parameter in the CO fits. This re- 
vealed that _R is a decreasing function of slit size for both 
the E140H and the E140M gratings of STIS (see Fig. 4 of 
iSheffer et al.ll200"7f ). The range of fitted resolving powers 
pe r aperture agr ees well with the range of values given 
bv lBoweri ()1997() . showing that Ismod.f has a good han- 
dle on R. Table 4 provides logarithmic values for total 
line of sight N of all five diatomic molecules that were 
modeled with Ismod.f spectrum synthesis in this study, 
as well as supplementary N values from the literature. 
Throughout the paper, all log N values are expressed in 
units of cm~^. 

3. OBSERVED RELATIONSHIPS 

3.1. Component Structures 

Our aim of deriving accurate column densities was the 
main reason for obtaining high-resolution spectra, both 
for CO from HST UV exposures and for CH, CH+, and 
CN from new optical data. These spectra reach veloc- 
ity resolutions between 1.5-2.0 km s~^, enough to re- 
solve many sight lines into multiple cloud components. 
The unveiling of such structures is important both for 
correctly treating the optical depth along the sight lines 
(by deriving b values for the line widths) and for distin- 
guishing characteristics of individual parcels of gas that 
would otherwise be lost in integrated line of sight values 
(|Pan et al.ll2005h . 

Our earlier anal ysis of indiv i dual c loud components to- 
ward Cep 0B2 in iPan et all (|2005[ ) showed that fits of 
CO cloud structures consistently resulted in structures 
that were similar to CN structures, even though the in- 
put for the synthesis of CO was based on CH cloud struc- 
tures that have more components than those found for 
CN. However, the sight lines sampled toward Cep 0B2 
were molecule rich, with a median value of A^(CO) = 
2.5 X 10^^ cm"^, and with CN detections along 73% of 
the sight lines. The new CO sample presented here is 
molecule poor, with a median value of iV(CO) = 1.0 x 
10^"* cm~^, i.e., a factor of 25 lower in CO abundance 
relative to the Cep 0B2 sample. Most of the poorest 
sight lines here (log A^(CO) < 14) are without detected 



A^(CN) but with CO cloud structures that are very simi- 
lar to those of both CH and CH+, whereas for sight lines 
with CN detections, CN is found in only about half of 
the components detected in CO. 

The lower molecular abundance along the new sight 
lines is also reflected in the generally low number of 
cloud components for all observed molecular species, 
even though these are mostly sight lines with large path- 
lengths that have a higher chance of intersecting molecu- 
lar clouds. Whereas 53% of the sample stars are farther 
than 1 kpc, and 23% are farther than 3 kpc, the derived 
cloud structures have low means of components per sight 
line: 1.9 ± 1.0, 1.8 ± 0.9, and 1.7 ± 0.8 for CH+, CH, 
and CO, respectively. This shows that CO is in excel- 
lent agreement with both CH and CH"*", underscoring the 
prevalence of low-density gas along these sight lines. The 
mean of CN components per sight line is smaller, 1.2 ± 
0.4. CN is detected along sight lines with the higher val- 
ues of iV(CO) and A^(H2), and then, on average, inside 
a single component. Among the 15 new CO sight lines 
with CN detections, 13 (87%) have N{GO) > 14.0, which 
is the median value for this sample. There are only three 
CO components that have CN with no detected CH+, 
and seven that are associated with CH+ but not with 
CN. We shall see in § 4.3 that this dichotomy between 
A^(CN) and Ar(CH+) can be employed as a good quali- 
tative indicator of nu- 

In the rest of this section we explore the correlations 
among logarithms of observed column densities by de- 
riving power-law parameters from regression analyses 
of the form log JV(My ) = log A + B x log A^(Mx) 
(jFederman et al.l Il990f ). where Mx and My are two 
molecular species. Unless otherwise indicate d, our BCES 
least-squares fits (jAkritas fc Bershadvlll996D are done on 
detections only, excluding the small number of upper lim- 
its. 

3.2. CO versus H2 

Panel (a) of Fig. 6 shows 105 sight lines with CO and 
H2 dete ctio ns taken from our sampl e and from t he sam- 
ples ofJCrennv& Fcdcr-mai^ (|200l . IPan et all 11005), 
and ISheffer et al.1 (120071). as well as inc l uding results to- 
ward bright stars from IFederman et al.l (|2003[ ). A single- 
slope global correlation returns a slope oi B = 1.89 ± 
0.15, having a correlation coefficient r = 0.834 and confi- 
dence level (CL) > 99.99%, see Table 5. The first indica- 
tion of a global correlation between iV (CO) and A^(H2) 
was provided bv IFederman et aTl ()1980l ). who found a B 
of «2 for 19 < log 7V(H2) < 21. Our fit of the new sam- 
ple also confirms the result of lLiszt fc Lucas! ()1998f) . who 
plotted N valu es of CO and H 2 from an updated version 
of the IFederman et al.l (| 19941 ) compilation to derive B = 
2.0 ± 0.3. 

From the start, the need for a variable slope descrip - 
tion of CO versus H2 was present. IFederman et al.l ()1980l ) 
found that lower- TV sight lines have a shallower B of »1.5, 
implying a slope break between 20.0 < lo g iV(H2) < 20.6 
and b etween 13.6 < log N{CO) < 14.8. iRachford et"an 
()2002l ) were the first to analyze higher-J V(H2) sight lines 
from FUSE, together with the data from lFederman et al.l 
(1994). They showed qualitatively that the slope of H2 
versus CO becomes shallower above log iV(H2) ~ 20.5. 
According to our inspection of their Fig. 3, the slope of 
the CO versus H2 relationship appears to get as high 
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as ft^S.S. Anoth er exploration of this relationship by 
iPan et al] ()2005f ) was based on a sample of FUSE sight 
lines toward the Cep 0B2 and Cep 0B3 associations. 
Despite gas density differences, these two associations 
presented similar CO/H2 slopes that also indicated a 
steeper relationship for log iV(H2) > 20, i.e., B = 3.2 
± 0.3 and 2.9 ± 0.6 toward Cep 0B2 and 0B3, respec- 
tively^ 

iBurgh et al.l (|2007t l examined 19 s ight lines and plotted 
these toget her with t he results of Crennv fc FedermanI 
(|2004s) and lPan et"^ (|2()05;) . confirming that CO versus 
H2 w as described by a relations hip with B k, 2. Like- 
wise, ISonnentrucker et "all ()2007| ) agreed that the over- 
all appearance of CO versus H2 is as steep as found by 
iFederman et al.l ()19 80') , but that it also appears to have a 
steeper increase of CO with H2 for A^(C O) > 10^^ cm~^ , 
in agreement wit h the fi ndings o f Rachfo rd et al.l (|2002f ) , 
iPan et all (|2005f ). and of iBurgh et al.. (,2007^ 

Overall, the indications are that the log iV(CO) ver- 
sus log A^(H2) relationship is not strictly linear (single- 
sloped) but that the slope itself (i.e., the exponent B) is 
also a function of Af(H2). Indeed, when we restrict the fit 
to the lower end of the distribution, the returned slope 
is shallower, while the higher end reveals a steeper slope. 
In Fig. 6(a) we also show the 10-point means of the CO- 
II2 sample, revealing a clear signature of two slopes, or 
a dual power-law correlation between CO and II2. Us- 
ing the two slopes to solve for their intersection point, we 
find the break between slopes occurs at log A'^(H2) = 20.4 
± 0.2 and log iV(CO) = 14.1 ± 0.1. Finally, employing 
the break location, we fit the two resulting sub-samples 
to find two highly significant (>4 a) correlations (Table 
5). Thus below the break in slope, log iV(CO) oc (1.46 
± 0.23) X log iV(H2), while above it, log iV(CO) ex (3.07 
± 0.73) X log Af(Il2). These two slopes are in excellent 
agreement with previous estimates, confirming all indi- 
cations that a steeper slope was needed for sight lines 
with higher values of N . 

The global behavior of iV(CO) versus iV(Il2) may be 
understood better when the UV data is complemented 
with N values for 293 dark clouds detected by millimeter- 
wave CO emission, an d taken from the compilation of 
IFederman et all ()1990l ). We note that iV(Il2) values for 
dark c louds were not direc tly observed, but were in- 
ferred (jFederman et al.lll990D from the corresponding vi- 
sual extinction {Ay), which is due solely to II2 in these 
clouds. This procedure cannot be applied to diffuse 
clouds (where Ay < 5 mag) because the (hydrogen) 
molecular fraction, /(H2) = 2iV(H2)/[A^(H I) + 2iV(H2)], 
is <1. As seen in panel (b) of Fig. 6, beyond the high- 
est end of the diffuse molecular cloud distribution one 
encounters the dark clouds, which have higher values of 
molecular N and of total gas density. 

The previous version of this connection between 
cloud cl asses was based on on ly 20 diffuse cloud data 
points (IFederman et al.l |1990D . At the time of the 
IFederman et al.l (|1990") study, only five CO values were 
known above log N = 15, with none available above log N 
— 16. At this time, however, we have 22 data points with 
log N{CO) > 15.0, of which five are above log 7V(C0) 
= 16, and one (HD 200775) is higher than log TV = 17. 
Thus during the intervening 18 years the gap between dif- 
fuse sight lines and dark clouds has been filling up with 
observations. [Simultaneously, recent radio observations 



are becoming more sensitive in their ability to measure 
smaller CO column densities that a p proac h the diffuse 
cloud regime, e.g., iGoldsmith et al.l (|2008D .] Not only 
are the two distributions seen to be stretching toward 
each other, but the steeper slope of the CO versus II2 
distribution, when extended to higher iV, is seen to pass 
near the center of the dark cloud distribution. However, 
the rise of -/V(CO) versus N(H2) cannot increase with- 
out limit because the supply of atomic carbon for CO 
formation will be exhausted. Thus the highest possible 
CO/H2 ratio is set by C/H2 = 2 x C/H or (2.8 ± 0.4) 
X 10~^, based on t he gas-phase C/H abundance from 
ICardelli et al.l (|1996f ) and shown in Fig. 6(b) as a dou- 
ble dashed line enclosing the ±1 ct range. Only a single 
datum out of 398 is seen to be slightly above the C/H2 
limit. 

One may need to go beyond fits with linear logarithmic 
slopes in order to allow for a better description of the CO- 
H2 relationship. Our test of a fourth-order polynomial 
fit returned a continuously variable slope that ranged 
from B = 0.8 to 3.4 for the lowest to highest iV(CO) 
values in the sample of diffuse clouds, at which point the 
slope declined while "connecting" with the dark cloud 
distribution. Such higher-order fits can only approxi- 
mate the more realistic slopes predic ted by d etailed CO 
photochemistry models. Bo th Rachf ord et al.l (2002) and 
ISonnentrucker et al.l ()2007f ) noted the agreement between 
the observed trend of CO versus H2 and results from the 
CO p hotochemical modeling of Ivan Dishoeck fc Blacld 
(jMl) for the CO-rich sight lines with log A^(CO) > 15. 
In Fig. 7 we show that the theoretical models of "translu- 
cent clouds" from Ivan Dishoeck fc BlacM ()1988f ) provide 
a functional variation that closely mimics the observed 
distribution of both diffuse and dark clouds, as well as 
along the transition region (i.e., their relative locations 
on the plot). These curves have steeper slopes than those 
obtained in our fourth-order polynomial fit, reaching as 
high as B = 4.3, 5.4, or 7.3 for the /uv = 0.5, 1.0, or 
10 models, respectively, where /uv denotes the enhance- 
ment factor over the mean interstellar UV radiation field. 
All three curves end up at the highest column densities 
with B between 1.5-1.9, i.e., bracketing the dark-cloud 
slope of i? = 1.62 ± 0.07. Our modeling with Cloudy 
of the CO versus H2 relationship that results from CH+ 
chemistry will be described in § 6. 

3.3. CH versus H2 (and CO versus CH) 

Based on 19 data points, iFedermaiil (|1982l ) demon- 
strated that A^(CH) is p roportional t o JV(H 2)i finding 
a slope of 1.0 ± 0.1, while iDanks et all (|1984[ ) confirmed 



this result by finding a slope of 0.85 ± 0.15 based on 
a slightly larger sample with lower-S/N data. (A com- 
bination of the two sa mples resulted i n a sl ope of 0.90 
± 0.10.) As for CO, lRachfordet"an (|2002D presented 
this relationship qualitatively, confirming its nearly lin- 
ear appearance and it s agre ement with the models of 
Ivan Dishoeck fc Blackl ()1989[ ) for the highest values of 
A^( CH) and N(R2)- A slope of 0.95 ± 0.10 was found 
bv IPan et all ()2005l ) toward 11 stars in Cep 0B2, but 
the only four data points that were available from the 
Cep 0B3 sample did not provide a clear case for a CH 
versus H2 correlation. Our log-log plot of CH relative 
to H2 (Fig. 8) shows a wcll-corrclatcd 90-point sample 
with a single slope oi B ^ 0.97 ± 0.07. Thus CH is defi- 
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nitely linearly related to H2, but as is the case with CO, 
the width of the correlation is appreciably larger than in- 
dividual measurement uncertainties. These correlations 
have CL above 99.9%, and thus our methodology used 
above in § 2.5, of importing CH cloud structures into 
spectrum syntheses of H2, is vindicated. 

Consequently, the ratio CH/H2 is a quantity that 
shows no correlation with H2. Our sample average is 
CH/H2 = 3.5;2i X 10-^ or log CH/H2 = -7.46 ± 
0.21. This value ca n be seen to agree well with the data 
plotted in Fig. 2 of lFedermaiil (| 19821) . Thus iV(H2) can 
be predicted from optical observations of A^(CH), pos- 
sibly with the exception of certain prominent photon- 
dominated, or photo-dissociation, regions (PDRs). The 
PDR targets HD 34078 and HD 37903 are found on the 
outskirts of the distribution, deviating from the average 
by about -1-3 and —2.5 u, respectively. 

It would be interesting to include CH and H2 val- 
ues for dark clouds to see how they rel ate to the plot - 
ted distribution of diffuse sight lines. iMatt ila' (1986) 
found that A^(CH)/7V(H2) = lO"'^'' for a smaU sample 
of dark clouds, i.e., in excellent agreeme nt with our av - 
erage above for diffuse molecular clouds. iMattilal ()1986l ) 
compared CH versus H2 for both types of clouds and 
found that they are in complete agreement, the dark 
cloud data being a monotonic extension of the CH-H2 
diffuse cloud relatioship, with a global slope oi B = 1.02 
lb 0.04. This i s confirmed here with the inclusion of the 
IMattilal ()1986l ) sample of dark clouds into Fig. 8, yield- 
ing B — 1.03 ± 0.03. [This relationship breaks down, 
however, for den se molecular c loud cores with N{]l2) ^ 
3 X 10^2 cm-2 ()Mattilalfl986[) .] In summary, the hn- 
ear relationship between CH and H2 for diffuse and dark 
molecular clouds is bound to be very useful for determi- 
nations of A^(H2) along sight lines where no CO data are 
available, or to corroborate such determinations based 
on CO data, as we shall show in § 4.4. 

The tigh t correspondence be t ween CH and H2 was em- 
ployed by iMagnani &: Onellol (|1995D in order to derive 
A^(H2) from millimeter-wave detections of A^(CH) toward 
diffuse and dark molecular clouds. Thus they were able 
to find variation by a factor of 20 in Xqo for diffuse 
clouds {Av < 4). Another relationship between CH and 
E(B — V) was presented bv IMagnani et al.l ()2003[ ). show- 
ing that while A^(H2) can be predicted from millimeter- 
wave CH and reddening measurements, CO cannot be 
derived from linear relationships with these parameters. 
Reddening values were already known to correlate with 
the total prot on density N(il) = 7V(H I) + 2A^(H2) in 
diffuse clouds fBohlin et al."1978! ), and with i V(H2) along 
translucent sight lines (Rachford et al.ir2002f ). 

The linear relationship between CH and H2 means 
that a plot of CO versus CH should be similar to 
the plot of CO versus H2. Indeed, our single-slope 
fit of 92 data points shows that B = 2.05 ± 0.21, 
i.e., A^(CO) varies as the sq uare of 7V(CH). Originally, 
iFederman fc LambertI ([1988) fitted a sample of 19 data 
points to find B = 1.97 for CO versus CH. Employing 
a sample twic e as large, this quadrat ic relationship was 
confirmed bv IFederman et al.l (|1994f ). who commented 
that above about iV(CH) ~ 3 x lO^^ and iV(CO) ~ lO^^ 
cm~^, A^(CO) is increasing more rapidly (higher B). In 
our sample, which is more than double in size yet again, 
one can see (Fig. 9) that the slope appears to be in- 



creasing above log A^(CH) - 13.3 and log iV(CO) - 15 
towa rd the locus of HD 200 775 in very good agreement 
with IFederman et al .! fl994') , who also noted the outly- 
ing position of the PDR sight line toward AE Aur (HD 

3407 8). 

Sonnentrucker et "all ()2007[ ) found a single slope of _B = 
4.0 ± 0.3 for a CO versus CH sample that tended to have 
higher column densities, i.e., mostly with log A^(CO) > 
14 and log A^(CH) > 13. In fact, as shown in Fig. 9, 
the 10-point means of our sample reveal the presence of 
two slopes with CL > 99.99%, B = 1.50 ± 0.30 and 2.80 
± 0.85, below and above a break at log iV(CH, CO) = 
(13.0, 14.1), respectively. There is excellent agreement 
with the power-law break found for CO versus H2 , since 
according to the CH/H2 ratio found above, log A^(CH) 
= 13.0 corresponds to log 7V(H2) = 20.5. As is the case 
for CO versus H2 , the steeper power law can be extended 
to reach near the center of the dark cloud distribution, 
which was plotted after converting its H2 coordinate into 
a CH location. These characteristics of the CO versus 
CH plot confirm that CH can be used as a dependable 
proxy for H2, and that a distinct change in CO photo- 
chemistry occurs at log A^(CO) = 14.1 ± 0.1. 

3.4. CH~^ versus H2 (and CH^ versus CO) 

iFedermaiil ()1982D also compiled and plotted 25 sight 
lines with detected CH+ and log A^(H2) > 19, finding an 
insignificant co rrelation (r = . 3 and CL <90%) with an 
unspecified B. iRachford et al.l (|2002[ ) also found a linear 
relationship between the two species, but with much in- 
crea sed scatter ab ove log A^(H2) ^ 20. They commented 
that lGredell (fl997l) found iV(CH+) and A^(CH) to be cor- 
related, which in the light of the tight correlation between 
CH and H2 presented in the previous section is consis- 
tent with a correlation between CH+ and H2. Indeed, 
our sample of 86 points returns a single-slope B — 0.42 
± 0.10 with CL > 99.99%, confirming that a global cor- 
relation exists between CH+ and H2. 

However, in agreement with lRachford et al.l (|2002f) . the 
current sample also shows a marked increase in data scat- 
ter or a loss of correlation above iV(H2) « 2. x 10^° 
cm~^. Panel (a) of Fig. 10 shows that employing 10- 
point means reveals yet another broken-slope relation- 
ship, with B = 0.78 ± 0.22 for log iV(H2) < 20.3, and 
B — 0.15 ± 0.21 above that break. The steeper power 
law indicates (at CL = 99.95% or 3.5 a) that CH+ varies 
nearly linearly with H2 at lower column densities. Above 
the break in slope, which is in excellent agreement with 
the H2 break from Fig. 6, both B and its CL indicate 
t hat CH+ and H2 ar e no lo nger correlated. 

iLambert fc Danka (jl986l ) presented correlations be- 
tween log A^(CH+) and log A^(H2*), i.e., column densities 
of excited states of H2 involving the J = 3 and 5 levels. 
Here we are unable to confirm these findings, because for 
our sight lines we do not find any CH+ and H2* corre- 
lations fo r all J = 1 to 4 lev els. However, we note that 
the ILambert fc Danka ()1986t ) study included sight lines 
with a range of N values much larger than ours, includ- 
ing sight lines that are H2-poorer by at least 3 orders of 
magnitude than those studied here. We shall comment 
further about excited H2 in § 6.2 when discussing the 
formation of CH+ . 

When CH+ is plotted against CO, the behavior is sim- 
ilar to that found just above for CH+ versus H2, namely. 
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presenting a slope break that is flanked by two different 
power-law fits (Fig. 10b). Below log N{CO) = 14.1 we 
find B = 0.46 ± 0.10 (CL > 99.99%) while above the 
break B = -0.14 ± 0.07 (CL = 93%). Again, showing 
excellent agreement with the CO break in Fig. 6, CH+ 
is definitely correlated with CO for lower column densi- 
ties, but is insignificantly (1.8 cr) anti-correlated with CO 
above the break. We believe that this behavior reflects 
the importance of CH+ reactions when 7V(C0) is low, 
so that a different production route operates at higher 
A^(CO). It is interesting to note further that with log 
7V(CH+) ex Bi X log iV(H2) and with log iV(CH+) ex B2 
X log A^(CO), one expects log A^(CO) to be ex B1/B2 
X log N{H2). Thus 0.78/0.46 = 1.7 ± 0.6 is in excel- 
lent agreement with B = 1.46 ± 0.23 found in § 3.2 for 
the low-A^ sight lines. Finally, both CH+ breaks relative 
to H2 (Fig. 10a) and to CO (Fig. 10b) are in complete 
mutual agreement that the abundance of CH+ presents a 
power-law change at log N = 13.2 ±0.1. A more detailed 
and chemically-motivated treatment of CH+ and its re- 
lationships with H2 and CO will be provided in § 6 based 
on numerical modeling with Cloudy and the incorpora- 
tion of a nonequilibrium term in the chemical formation 
of CH+. 

3.5. CN versus H2 (and CO versus CN) 

CN detections here encompass a smaller sample that 
includes less than half of the sight lines t hat are included 
in the H2, CO, CH, and CH+ samples. iFederman et al.l 
()1984[ ) analyzed a smaller sample still, and concluded 
that the CN abundance wa s prop ortional to the third 
power of H2 . iRachford et al.l (|2002f ) showed that a strong 
correlation exists here as well, with an estimated B ^ 2.5 
according to our inspection. Such steep slopes are not 
confirmed here, because we find B = 1.5 ± 0.4 from our 
regression fit, with r = 0.67 (CL > 99.99%). However, 
owing to the relatively large number of upper limits on 
CN, we decided to employ the Buckley- James method of 
linear regression with censo red data, available from the 
ASURV statistical package (jlsobe et al.lll986[) . This fit 
returned a steeper slope of i3 = 1.8 ± 0.4 (Fig. 11), 
which was also confirmed by the EM algorithm and by 
Schmitt's method, but is also appreciably shallower than 
B K, i. One possible explanation for the disagreement 
among fitted slopes may involve the smaller number of 
sight lines with detected CN in previous studies. Our 
sample includes 40 sight lines with detected CN, and as 
a result there is a significant "re-" population of the plot 
with log Af(H2) < 20.5. However, the sample is still 
too small and restricted in range to reveal information 
about any power-law break. Thus arbitrarily breaking 
the sample at mid-point results in similar slopes, both of 
which have CL below 90% (Table 5). 

Figure 12 shows the run of log A^(CO) versus log 
Af(CN) with B = IA± 0.2, fitting the CN detections 
only. This 42-point sample has r = 0.84 and thus CL > 
99.99%. Again, owing to the presence of 17 CN upper 
limits, we employed the Buckley- James censored-data 
fit and derived B = 1.9 ± 0.2, a slope that seems to 
treat the upper limits as quasi-detections. On the other 
hand, using the Schmitt method returned a fit with the 
significantly lower B = 0.8 ± 0.2, which seems to ex- 
clude all upper limits in the independent variable (i.e., 
in CN) . Thus our result based on detections only is brack- 



eted by the two slopes based on censored-data methods. 
A smalle r sample of (uncensored) 31 points was pre- 
sented bv lSonnentrucker et al.l (|2007D . yielding (in an un- 
weighted fit) S = 1.5 ± 0.1, and r = 0.90 (CL > 99.9%), 
thus agreeing more with our detections-only fit than with 
our Buckley- James fit. Again, splitting the sample into 
two equal sub-samples returns two <2 a fits with iden- 
tical slopes, i.e., with no evidence for a break. Since the 
CN sample is derived for the most part from high-iV sight 
lines that are above the CO versus H2 break, consistency 
is still preserved. 

4. DERIVED PHYSICAL CONDITIONS 

4.1. UV Shielding of CO 

Earlier (§ 3.2) we described the finding of a power- 
law break in the correlation analysis of CO versus H2 
at log N — 14.1 versus 20.4. This break in slope 
was also confirmed above through analysis of CO ver- 
sus CH, CH+ versus H2, and CH+ versus CO. This 
power-law break corresponds to a change in CO pho- 
tochemistry. Our value is sim ilar to log A^(CO) > 14 
found by Frcrk ing et al.l (|1982f ) for two CO isotopologues, 
C^^O and ^'^CO, through their correlations with Ay in 
molecular clouds. In addition, a value of log 7V(C0) 
> 14 was supported by t he self-shielding computations 
of iBallv &: Langen ()1982f ). who showed the importance 
of line photodissociation in steepening the increase of 
A ^(CO) with depth into a mole cular cloud. 

Ivan Dishoeck fc Blackl ()1988f ) presented the most de- 



tailed modeling of CO photochemistry in the regime 
of translucent (1 < Ay < 5 mag) sight lines. In § 
3.2 (and Fig. 7) we also compared the observed dis- 
tri bution of diffuse and dark cl ouds with three families 
of Ivan Dishoeck fc Blackl ()1988l ) models for translucent 
sight lines differing in /uv- This showed the good global 
agreement of ab undance trends for CO versus H2. In 
their modeling, Ivan Dishoeck fc Blackl (jTgSS) incorpo- 
rated a detailed description of numerous CO absorption 
bands in the far-UV, owing to their importance in dimin- 
ishing CO photodissociation rates through self shielding. 
This detailed band structure was also needed for precise 
accounting of the shielding of CO bands by H2 absorp- 
tion lines, since the total UV shielding of CO is con- 
trolled by both A^(CO) and A^(H2). This two-parameter 
shielding funct ion of CO , 8, w as tabulated in Table 5 
of van Dishoe ck fc Blackl (|1988D . The values of 6 are 
smaller than 1 since they provide the reduction in the 
photodissociation rate of CO. When O < 0.1, or log 
A^(CO) > 15, total UV shielding of CO rcsuffs in rapid 
steepening of the CO versus H2 relationship. 

In order to determine the effect that 8 has on the ob- 
served distribution of CO versus H2, we present again 
in Fig. 13 the global view of the diffuse and dark cloud 
distribution, together with our dual power-law fits, all 
overlaid by contours of th eoretical 8 v alues based on an 
interpolation of Table 5 of lvan Dishoe ck fc Black (1988). 
It can be seen that for decreasing values of 8 the values 
of A^(CO) are increasing, as the photodestruction of the 
molecule is being diminished. T his provides a demon- 
stration of the applicability of the Ivan Dishoeck fc BlacM 
(JT988) shielding function to diffuse and dark sight lines. 
Quantitatively, the location of the slope break at log 
A^(H2, CO) = (20.4, 14.1) is seen to be near 8 « 
0.4, i.e., where the reduction in the photodissociation 
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of CO equals 1 mag. Here, we interpret the slope break 
as the locus of a transition between lower-density and 
higher-density regimes, but the comparison with O sug- 
gests some contribution from the UV shielding of CO 
to the steepening of the slope near the break. Fu- 
ture calculations of O, which should include updated 
/-values of predissociating CO bands (iFederman et all 
[200lt iSheffer et al.l [20031 : lEidelsberg et al.l 12004 I2006D . 
may clarify the association between O and the observed 
break in the slope of CO versus H2 . 

4.2. Excitation Temperatures 

Each molecule in the ISM is influenced by both colli- 
sions with other molecules and atoms (matter) and by 
interactions with photons (radiation). The former pro- 
cesses will tend to thermalize the internal level popula- 
tions of the molecules so that their Tex will reflect the 
kinetic energy of the colliding gas particles. Such a case 
is reflected in the J = and 1 populations of H2 ow- 
ing to the lack of permitted dipole transitions that lead 
to cooling of the molecule. As can be seen in Table 
6, H2 along all sight lines has relatively high values of 
Tc^{J = 0,1) (hereafter Tqi) since they reflect the pre- 
vailing kinetic temperature of the gas. The average of 56 
sight lines with newly-derived H2 parameters is Toi(H2) 
= 76 ± 14 K, whic h is in excellent agreement with the 
iSavage et all (Il977f ) resul t of 77 ± 17 K. As we pointed 
out in ISheffer et aLl (|2007f) . sight lines with lower Tqi val- 
ues for H2 are associated with detected amounts of ^■^CO. 
The sub-sample of 25 sight lines with ^'^CO was shown 
to have an average Toi(H2) that is 20 K below the aver- 
age Tqi for sight lines without detected ^^CO. The entire 
sample here is not constrained by the presence of ^'^CO, 
just like the original sample of ISavage et al.l (119771 ). 

The YS Ismod.f code also returned fitted Tox values 
for the higher- J levels of H2, as listed in Table 6. A cur- 
sory inspection shows that for H2, as detected in diffuse 
molecular gas, Tqi < T02 < T03 < T04. This is confirmed 
by their means from all sight lines: 76 ± 15, 101 ± 15, 
140 ± 23, and 213 ± 31 K, respectively. As can be seen 
in Fig. 14, logarithmic correlations exist between each 
Toj (for J > 1) and Tqi, all with CL > 99.99%. The re- 
spective slopes of the regressions are mutually identical 
within their uncertainties: B = 0.48 ± 0.07, 0.52 ± 0.09, 
and 0.46 ± 0.08, in order of increasing J. The positive 
slopes may be indicating a connection between increas- 
ing Tqj and decreasing gas density at cloud edges, where 
there is more efficient pumping by FUV photons. 

CO presents the opposite case of sub-thermal excita- 
tion, with Tex values rarely rising above «5 K. The only 
interesting case of warmer CO is found along the p Oph 
D (HD 147888) sight line: Toi(CO) = 13.6 K. In this case 
the CO being probed is near the p Oph molecular cloud, 
and CO emission from the latter is able to raise Tex(CO ) 
along the diffuse part of the cloud (jWannier et al.llT997f ). 
The average Tqi (CO) from 61 (62) sight lines without 
(with) p Oph D is 3.5 ± 0.7 (3.6 ± 1.5) K. Slightly higher 
values are found for To2(CO) = 4.2 ± 0.8 (4.4 ± 1.4) K 
and To3(CO) = 5.3 ± 1.3 (6.0 ± 2.0) K, without (with) p 
Oph D, but more likely showing that Te2:(C0) is constant 
within the uncertainties. We also see no dependence of 
Toi(CO) on the density indicator CN/CH+. Further- 
more, the same mean value for Tqi (CO) is found for sight 
lines with or without detected CN, despite a density dif- 



ference of a factor of 10. These differences from Te2;(Il2) 
arise because in diffuse molecular clouds densities remain 
below the critical density of ~2000 cm~^ for CO. 

4.3. Total Gas Density: Empirical Indicators 

The total gas (proton) density, nn = '^(H I) -I- 2n(H2), 
controls the chemical reaction networks via the density 
dependence of molecule-production terms, and thus af- 
fects the resultant molecular abundances. We may ex- 
plore such density effects by deriving njj for sight lines 
in a variety of ways. The simplest and most empiri- 
cal invo lves the ratio of two observables, A^(CII) and 
N{CN). iCardein et all (Il99lh showed that the CN/CH 
ratio is correlated with n^ since CN is formed inside 
denser and colder clumps of gas out of pre-existing CH. 
Thus plotting other quantities versus CN/CH is tanta- 
mount t o showing the relationship of those quantities 
with nn (jSonnentrucker et al.ll2007D . 

As will be reinforced below, the CH+ molecule is typ- 
ically formed in lower density regimes, leading to a de- 
pendence 0££Osite_to_thatofCN. This was shown empir- 
ically by iCardelli et aLJ (|l99l who found that CN/CH 
was anti-correlated with CH+/CH. Thus, instead of us- 
ing CH, which is less dependent on densi ty thanks to its 
conn e ction to gas cont aining CN or CH+ (iLambert et al.l 
I1990I iPanet alilMol . CN/CH+ should be a more ef- 
fective indicator of the density than the CN/CH ratio. 
In Fig. 15 we show that both empirical density indica- 
tors are well correlated with each other (r = 0.77, CL > 
99.99%). Note, however, that CN/CH extends over less 
than two orders of magnitude, whereas CN/CH+ spans 
three orders, suggesting that CN/CH+ responds better 
to changes in rin, a picture consistent with the presence 
of CH in both high- and low-nn gas. 

Panel (a) of Fig. 16 shows the abundance of CO rela- 
tive to H2 versus the empirical density indicator CN/CH. 
With r = 0.60 the plot shows a very good correlation (CL 
> 99.99%) between the two quantities, having B = 1.16 
± 0.25. For comparison, Fig. 16(b) shows CO/H2 ver- 
sus CN/CH+, and the correlation is found to be even 
tighter, having a larger r of 0.78 and B = 0.85 ± 0.11. 
Similar plots of species other than CO also exhibit larger 
r values and visibly tighter relationships, confirming the 
better role of CN/CH+ in sorting diffuse sight lines ac- 
cording to rifj- Furthermore, the tighter correlation seen 
in Fig. 16(b) is an indication that the CO/H2 ratio is 
controlled significantly by the local gas density. 

Given the better association between nn and the ob- 
served CN/CH+ ratio, one may imitate a 3-dimensional 
(3-D) plot by employing proportionately-sized symbols 
to represent values of the latter quantity on the 2-D sur- 
face of, e.g., the A^(CO) versus iV(H2) plot that was 
shown in Fig. 6. From the resulting Fig. 17 one can 
discern two general trends involving variations in gas 
density (as given by CN/CH+). First, density is clearly 
the lowest toward the lower left corner of the plot (where 
many of the values are upper limits) and vice versa, show- 
ing that both A^(CO) and A^(H2) are correlated with nn- 
Second, the density clearly varies in an orthogonal di- 
rection to its first gradient, i.e., it is increasing from the 
lo wer envelope to the upp er envelope of the distribution. 

Fed erman et al.l ()1980( ) were the first to find that the 
dispersion in the relationship of CO versus H2 is signif- 
icantly larger than the measurement uncertainties asso- 
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ciated with individual data points, as confirmed in Fig. 
6. These authors were also able to show that applying 
high-density and low-density chemical models to this re- 
lationship indicated that its width was affected by gas 
density, such that nn is higher at the upper envelope, 
in agreement with our findings using CN/CH+ as the 
density indicator. 

The two Cepheus samples in iPan et all (|2005D ap- 
peared to occupy nonoverlapping positions on the CO 
versus H2 plot. However, when compared with the cur- 
rent, much_Jajgei_sample that includes the sight lines 
from lPan et al.l (j2005f ). all Cep 0B2 (having higher den- 
sity gas) and 0B3 (lower density) data points are part 
of the global distribution of points, although they seem 
to belong to the upper and lower envelopes of the distri- 
bution, respectively, thus providing more support to the 
overall pict ure. 

Recentlv. iLisztJ (|2007( ) suggested that the CO versus H2 
relationship in diffuse clouds directly reflects the forma- 
tion of CO from a HCO+ precursor. However, while the 
second (cross- wise) variat ion agrees qualitatively with 
the models of iLisztl ()2007() . see his Fig. 1, those mod- 
els do not reproduce the first variation, i.e., the rise in 
den sity in tand em with increasing N values. In fact, 
the lLisztl ()2007f ) models have constant density values be- 
tween the lower left and upper right corners of the CO 
vers us H2 plot. Perhaps this difference is an indication 
that lLisztl 's assumption that CO production is controlled 
by recombination of HCO+ with a constant abundance 
of 2 X 10~^ relative to H2 is inadequate. 

4.4. Predicted Column Densities for H2 

After fitting a correlation plot between two observed 
column densities (-/Vq), fit parameters may be used to 
predict one of the column densities (iVp) in the absence 
of the other. The value of ^p(H2) = 7.4-10.0 x 10^ ° 
cm-2 toward HD 208266 was given in lPanet all (|2005f ). 
based on their fits of ^o(CO) and iVo(CH) versus No(ii2) 
for the small sample of sight lines toward Cep 0B2. Here 
we use the dual-slope relationship between CO and H2 in 
Fig. 6(a), as well as the single-slope relationship between 
CH and H2 in Fig. 8, to predict (the unobserved) A^(H2) 
toward the two stars without H2 data but with CO and 
CH data, HD 36841 and HD 43818/11 Gem, which are 
near the CO versus H2 break in slope, as well as toward 
HD 208266. These predictions employ the global ±20% 
1-CT uncertainties in N^ values. 

The CO-based log iVp(H2) values for HD 36841 and HD 
43818 are 20.46 ± 0.06 and 20.32 ± 0.06, respectively 
The same exercise for HD 208266 yields log A^p(H2) = 
21.12 ± 0.03, which is 2 a away from the predicted range 
(20.87-21.00) given in Pan et al.l (|2005l. The CH-based 
log A^p(H2) values for HD 36841 and HD 43818 are 20.41 
± 0.08 and 20.45 ± 0.08, respectively, which are 0.05 
lower and 0.13 higher than the CO-based log Np{'il2)- 
Conceivably, when the difference is larger than 20% (or 
>0.08 in the log) it is reflecting the additional uncertain- 
ties introduced by the intrinsic widths of the correlations. 
The same exercise for HD 208266 yields log Nplji2) = 
20.95 ± 0.08, which is 0.17 lower than the CO-based 
prediction but in excellen t agreemen t with the predicted 
mid-range value given in iPan et al. (2005). Combining 
results from CO and CH, both HD 36841 and HD 43818 
are predicted here to have log A'p(H2) = 20.4 ± 0.1, while 



for HD 208266 the prediction is 21.0 ± 0.1. The corre- 
sponding 3-(T uncertainties are provided by the full width 
(±0.3) of the horizontal spread of iVo(H2) in both Fig. 6 
and Fig. 8. 

5. ANALYTICAL CHEMISTRY 

We examined the results presented here from two 
chemical perspectives. First, two sets of analytical ex- 
pressions from previous work are used in this section to 
extract njj associated with the material containing CO, 
one set involves the equilibrium chemistry leading to CN 
and another set describes the (equilibrium) synthesis of 
CH from CH^ . In § 6 we provide a more general chemical 
analysis based on the use of the Cloudy code. 

5.1. CN Chemistry 



Analytical expressions 
ing CH, C2 , and CN 



for the chemistry connect- 
diffuse interstellar clouds 



Federman et al.l I1994D. with update d rate coefficients 



upgate d 
Knauth et al.l 1200 It iPan et all 1200 H ). are used to ex- 



tract estimates for gas density, nH(CN). To sum- 
marize, the production of CN is primarily given 
by the reactions C2(N,C)CN, CH(N,H)CN, and the 
chain C+(NH,H)CN+(H2,H)HCN+(e,H)CN, which has 
a parallel path where HCN+ reacts with H2 to pro- 
duce H2CN+ and then CN via electron recombina- 
tion. Observed A^(CH) and iV(C2) (when available) are 
adopted for the comparison between observed and pre- 
dicted CN column densities. A steady-state rate equa- 
tion involving terms for chemical production and photo- 
dest ruction of CN is emp loyed in the determination of 
riH (jFederman et al.lll994D . 

As in our r ecent papers (e.g ., iGredel et al.l 120021 : 
IPan et al.l[2005l : IWeltv et al.ll2006l ). results are presented 
for i ndividual vel ocity components, whenever possible. 
Like iGredel et aLl (|2002), we determine upper limits on 
nH(CN) for components without detectable amounts of 
CN absorption. This is especially important for the re- 
sults presented here because many new detections of CO 
are found in directions that are not very rich in molecules. 
We do not repeat the chemical analysis for sight lines 
in Ophiuchus an d in Cep 0B2 and 0B3 described in 
iPan et all (|2005D . nor the analyses of the photodissoci- 
ation reg ions illuminated by H D 37903 and HD 200775 
found in ' Knauth et al.l (|2001l ) . Similarly, many of the 
directions contained in a rea nalysis of spectra acquired 
with the Coperni cus satellite (ICrennv fc Federmanll2004f) 
are discussed in IZsargo fc FedermanI ()2003[ K Updates 
are given for some of th e sight lines e xamin ed by 
iFederman et all ()1994D and IWannier et~all ()1999[ ) in or- 
der to provide a self-consistent analysis and comparisons 
are presented below. 

A key ingredient in this analysis is the value for 
the amount of extinction at UV wavelengths caused 
by interstellar grains, ryvj for each sight line. This 
was determined by examining various measures 
for grain pr operties: the ratio of tot al to selective 
extin cti on (|Cardelli fc Clavtonl 119911; iL arson et al.l 
1 20001: iBa^baro et al.l 120011: iP atriarchi et all l2001t 



'Whittet et al."2001': 'Ducati et aL 2003; Patriarchi et alj 
2003; Valcncic et al. 2004; Fitzpatrick fc Massa 200i^ 
Tarson fc Whitte t 2005; So fia et al.ll2005D. the shape of 
the UV extinction curve (.Massa et al.l 19831: IWitt et al.l 
119841: iFitzpatrick fc Massal 119901 : iPapai et all Il991l : 
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Weltv k Fowled [l99l iLarson. Whittet. fc Houehl \T9M. 



Patriarchi fc Perinottol Il999( ). an d a comparison of 
the r atio ^(15 - V)/E(B-V) (jK relowski fc Strobell 
[1981 iSavaee et alJflQSa iPaoai fc fc elowski 1992). For 
most sight Unes, typical grain properties apply and we 
adopted ryv = 2 x 3.1E{B — V), where 2 is a prefactor 
that depends on chara cteristics of the extinction curve 
(jFederman et al.1 I1994D and the amount of reddening 
came fr om the work cited above o r from ISeab fc Snow! 
(fT98l . iCarnochanI ((19861) . and lAiello et al.1 (|1988D . 
When grain properties suggested enhanced UV ex- 
tinction, we used prefactors of 3 (for the directions 
toward HD 12323, HD 15137, HD 36841, HD 163758, 
HD 185418, HD 198781, and HD 210121) or 2.5 (for 
HD 14434) instead of 2, depending on the severity of 
the difference from typical values. Several directions 
(HD 96675, HD 99872, HD 102065, and HD 124314) 
indicated below typical UV extinction; here the pref- 
actor was set to 1.7. For HD 93840, an intermediate 
value seems to be appropriate, and a prefactor of 1.85 
was adopted. For sight lines without information on 
grain properties (HD 24190, HD 30122, HD 137595, 
HD 190918, HD 192035, and HD 192369), the typical 
relation was employed. One further constraint was 
considered: components having separations less than 20 
km s~^ were assumed to arise from nearby complexes 
where shadowing would be present and each component 
would experience the full amount of extinction as a 
result. The lone exception was the direction toward 
HD 13745, where the components are 26 km s"'^ apart. 
The results for HD 30122 given in Table 7 show the 
effects that uncertainties in ruv have on nH(CN). In 
general, photodissociation is the dominant destruction 
pathway for the clouds in our study, and therefore 
uncertainties in ruv lead to inferred uncertainties of 
-30% in nH(CN). 

For many of the directions, N data for CH, CN, as 
well as CH+ used in the next section, were obtained 
as part of the present study. Much of the remaining 
da ta on A^(CH) an d JV(C N) come from the compilation 
of iFederman et al.l () 19941 ) , but there are a number of 
update s now available. F or ^ Per, we i ncluded the re- 
sults o f ICrane et aD ()1995[ ) for C H and of iLambert et al.l 
(fT995h for C2. The results of ILambert et all (|l995h 
for ( Oph were also used here. We adopted the re- 
sults of iKaczmarczvkl (|2000l) for the C2 colu r nn tow ard 
X Per. The CH results of JAndersson et all ()2002f ) for 
HD 99872, HD 115455, and HD13 7595 ar e inclu ded, 
as are the CH and CH+ resuhs of iGredell (|1997D for 
HD 114886. For the sight line toward HD 154368 we 
incorporated A^(CH) from W elty (2005, private com- 
munication) and iV(CN) fromlRoth. Mever. fc Hawking 
(fl9933 and iRoth fc Meval ([19951). For the gas to- 
ward HD 185418 and HP 192639, w e used the results 
from ISoim entrucke r et al.l (l2002l l2003l ) , supplemented by 

Since [Pan et"aLl (|2005[ ) 



those of [rhorburn et al.l 
did not consider directions without detectable amounts 
of CN in their analyses, we do so here for HD 208440, 
HD 208905, HD 209339, 19 Cep, and HD 217035A. For 
stars in Per 0B2 (40 Per, HD 23478, and HD 24190) 
we used our unpublished results. We also note that for 
gas toward o Per, X Per, and 62 Tau, component struc- 
ture is available for C N, but not C2. Since th ese species 
appear to coexist (e.g.. IFederman et al.lll994l ). we scaled 



the C2 results so that CN/C2 was the same for each com- 
ponent. Finally, as indicated in the Table, line-of-sight 
results are given for directions where component infor- 
mation is missing for C2 as well as CN. 

The results of this analysis appear in Table 7. For each 
cloud in a specific direction, we list the observed values 
A'^o(CH), iVo(C2), and A'^o(CN), and the predicted values 
-/Vp(C2) and -/Vp(CN) that best match the observations, 
the kinetic temperature (T), I\jy, ruv, and nH(CN). N 
values are given in units of 10"'^^ cm~^. Most calculations 
are based on /uv equaling 1 and on T = 65 K. The latter 
value is not critical because the results for uh are not very 
sensitive to T. For especially molecule-rich clouds and for 
some clouds studied by us in the past, lower values for T 
are adopted. The A^p values are generally in very good 
agreement, and are always within a factor of two, of the 
No values. 

5.2. CH+ Chemistry 

For many of the directions listed in Table 7, only upper 
limits on CN are available. For nearly all of these, CO 
production via reactions involving CH+ appears likely 
(see below). We therefore considered estimating the 
gas d ensity from the chemical scheme link ing CH and 
CH+ (I Weltv et al.ll2006t iRitchevet aP [20061 ^ as follows: 
CH+(H2,H)CH+(H2,H)CH+ and the dissociative recom- 
bination CHJ'(e,H2) CH. In particu lar, we used the ana- 
lytical expression in lRitchey et al.l . 



nn 



A^(CH) 2/uvr(CH) 
A^(CH+) 0.67A:/(H2) ' 



where r(CH) is the CH photodissociation rate [1.3 x 
10~^ exp(— ruv) s^^], k is the rate coefficient for the 
reaction CH+(H2,H)CH+ (1.2 x lO'^ cm^ s'^), and 
/(H2) is the molecular fraction. In ad dition to the 
present study, A^(H2) v alues come from ISav age et al.l 
(191^, iRachford et"al[ (|2002f l. and iPan etal. ' (2005A. 
The colur nn densities o f atom ic hydrogen, A^(H I), 
are from ISavage et al.l (I1977D fo r the brigh t stars, 
and from S lachfo rd et al.l (l2002f). lAndre et al.l ([2 0031 , 
ICartledge et al.l (|2004f ). and iJensen fc Snow! ([ 20071 for 
sight lines studied with FUSE. Data on atomic hydro- 
gen do not exist for the sight lines toward HD 114886 
and HD 137595. For these stars, we estimated 7V(H I) 
from E(B — V) using the relationshi p between reddenin; 
and total proton column density of Bohlin et al.l (|197' 
and accounting for the amount of H2. Values of N for 
the carbon-bearing molecules are taken from the sources 
given in the previo us section for the mo st part or from 
those compiled by ICrennv fc FedermanI (2004) for the 
bright stars. The results appear in Table 8. 

5.3. Comparison of Results 

Many of the sight lines listed in Table 7 were analyzed 
in our previous work. For Per, X Per, and 62 Tau, we 
now incorporated results for individual velocity compo- 
nents by scaling the values for N{C2) to those measured 
for A^(CN). The updated chemistry does not sign ificantly 
affect the conclusions of IFederman et all ([1994f ) for the 
three sight lines, nor for the gas toward C Per, £, Per, 
X Oph, and C Oph. There is also reasonable correspon- 
dence be tween the prese n t chem ical results for 40 Per and 
those of IWannier et al.] (|1999[ ). which are based on C2 



UV Survey of CO and H2 



11 



(and C I) excitation. Finally, our inferred density for the 
main component toward HD 1543 68 is about a factor of 2 
larg er than our previous estimate ([Federman fc LambertI 
|198§). which was based on an earlier, higher measure for 
7V(CH). This refined value for nnfCN) is consistent with 
the results of Ivan Dishoeck fc de Zeeuwl fi 1984^1 from C;:> 
excita tion (300-1000 cm~^) and iBlack fc van Dishoeckl 
(|1991l ) from CN excitation (750-2000 cm'^). 

Comparisons for several other sight lines with other 
work are also available. The most well-studied of these 
directions is toward HD 210121. Our chemical analy- 
sis indicates that nn = 1425 cm"'^, one of the highest 
values in Table 7. The large value agrees with results 
from analyse s of m olecular excitation. For C2 and CN, 
iGredel et ai] (|1992f) found densities of 500 to 1000 cm~^ 
and 1 500 to 2500 cm~'^, respectively, while iRoueff et all 
(|2002| ) obtained densities of about 2000 cm~^ from the 
distribution of C3 levels. Our upper limits on density 
for the gas toward HD 185418 and HD 192639 (< 30 to 
40 cm~'^) are consisten t with de nsities inferred from C I 
excitation ( Sonnentrucker et al.ll2002. 200 3). 

For most directions with CN upper limits, results 
from CN cheniistry a re not very constraining (e.g., 
iFederman et al.lll997bf ). Instead, comparisons with C I 
excitation seem to be more appropriate, as in the cases 
of HD 185418 and HD 192639. For such directions, CH+ 
chemistry usually dominates, and so the results from Ta- 
ble 8 would be more meaningful. For the nearby bright 
stars, a Ca m, v Sco, and m Nor, a comparison w ith the 
results from lJenkins. Jura, fc LoewensteinI (|1983l ) is pos- 
sible. For the three directions, densities are derived from 
the quoted pressures assuming that Tqi is the kinetic 
temperature. Consistency between results occurs for a 
Cam and n Nor, where upper limits from C I excitation 
are 10 and 370 cm~'^, respectively, versus our values of 
about 1 for no enhancement in the strength of the UV ra- 
diation field. However, the lower limit toward v Sco of 70 
cm~^ contrasts with our density of 10 cm^^. If the cloud 
were near the star, such that /uv was greater than 1, 
the two measures could be brought into agreement. This 
point s out a deficiency in the current analysis of CH+-like 
CH (Lam bert et al.|[l990 ) for many of the sight lines in 
Table 8: Most stars in the Table lie at least a kpc away. 
Thus, a significant amount of atomic hydrogen and a cor- 
responding amount of UV extinction are not likely asso- 
ciated with the gas containing CH and CH+. The most 
important factor is the exponential change in extinction, 
thereby increasing the estimate for nu- The typical in- 
crease is about a factor of 10, but can be significantly 
more. For the clouds toward HD 185418 and HD 192639, 
densities approachi ng 10 are possible, comparabl e to the 
results from C I (S onnentrucker e t al. 2002l l2003f ). Over- 
all, the analytical results presented in Tables 7 and 8 are 
consistent with the detailed models described next. In 
particular, the gas rich in CN and CO tends to have 
larger densities. 



5.4. The Role of C2 Data 

Solving for nH(CN) via the analytical chemistry 
scheme depends on input of A^(CH) and 7V(CN), but 
not necessarily on the availability of 7V(C2). When 
A^(C2) is not available, it is predicted by the model 
via CH(C+,H)Cj, which is followed by Cj(H2,H)C2H+ 



and C2H+(H2,H)C2HJ, both leading to the C2 molecule 
by dissociative recombination. The predicted value of 
N{C2) is then determined from a steady-state rate equa- 
tion involving terms for C2 chemical production and pho- 
todestruction. However, when A^(C2) is known, it intro- 
duces a constraint on the chemical formation of C2 , which 
in turn is affecting the production (and predicted abun- 
dance) of CN. It would be interesting to investigate the 
differences, if any, between sight lines with or without 
observed values of A^(C2). 

Figure 18(a) shows a logarithmic plot of CN-derived 
gas densities, nH(CN), versus observed 7V(CN) for sight 
lines with detected CN components, fitted with B = 0.31 
± 0.10 (solid line). (Data wer e gathered from the chem- 
ical analys es of this paper, of 'K nauth et al.l|2001l and of 
iPan et clLI 12005.) Whereas sight lines without observed 
values of A^(C2) (empty circles) show a loose correla- 
tion with r = 0.32 between nH(CN) and A^(CN), a much 
tighter correlation with r — 0.61 is seen in the distribu- 
tion of the filled circles, which denote sight lines with 
C2 measurements. Thus it is apparent that more robust 
derivations of nH(CN) require not only observed A^(CH) 
and iV(CN) values, but also observed N{C2) data. Simi- 
larly, Fig. 18(b) shows the run of nH(CN) versus the ra- 
tio CN/CH"*", which is the better empirical proxy for gas 
density (§ 4.3). Again it is obvious that sight lines with 
observed iV(C2) show a tighter correlation than those 
without C2 data. 

Finally, as a confirmation of the close affinity between 
resultant gas densities and input C2 observations, we 
plot in Fig. 19(a) these two quantities. It shows that 
nH(CN) predictions correlate well with N{C2) for sight 
fines with C2 detections (r = 0.63, CL > 99.5%). This 
relationship, with B = 0.48 ± 0.15, is in line with the 
expectation that CN and C2 molecules are formed in 
higher density (and colder) clumps of gas because of the 
correlation found above between nH(CN) and iV(CN). 
Figure 19(b) indeed shows that the correspondence be- 
tween the observables A''(CN) and N{C2) has a slope of 
0.97 ± 0.28 and r = 0.64, or CL > 99.7%. This result 
ag rees at the 2-cr level wi th the slope of 1.6 ± 0.2 found 
bv IFederman et all (| 19941 ) from a larger 33-point sample 
that sho wed a ti g hter c orrelation (r = 0.85 and CL > 
99.99%). iGredell ()2004l ) showed that a small sample of 
sight lines toward Cep 0B4 is probing a high-nn molec- 
ular cloud and providing tight correlations between all 
three molecules that are involved in equilibrium chem- 
istry: CH, C2, and CN. It is clear that CN chemistry is 
dependent on A''(C2) and that robust nH(CN) values can 
be derived from these two observables. 

6. NUMERICAL MODELS WITH CLOUDY 

6.1. Computational Details 

We performed a series of model calculations designed 
to cover a range of physical parameters characteristic of 
diffuse and molecular clouds. Our calculations used ver- 
sion C07.0 2 of the spectra l synth e sis code Cloud y, last de- 
scribed bylFerland et aD (Il998[ l l^an Hoof et al.. (.20041 
lAbel et all (|2005f ). IShaw et all (120051 ). and lR511ig et all 
(|2007l ) discuss in detail the Cloudy treatment of vari- 
ous physical processes impor tant in modeling at omic and 
molecular phases of the ISM. Rollig et al. (2007) compare 
the predictions made by various PDR codes and find ex- 
cellent agreement between Cloudy and the codes used in 
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iKaufman et~aI1 (|1999D . iBoger fc Sternberd (l2005f). and 
the Meudon PDR code (|Le Petit et al.ll2006f ). 

The geometry of our model is a plane-parallel slab il- 
luminated from both sides by far-UV radiation. This 
geometry is appropriate for diffuse environments bathed 
on all sides by the far-U V radiation field, and is iden- 
tical to that used in both Ivan Dishoe ck & Blackl ()1988D 
and lLe Petit et all (|2006D . iLe Petit et al.i (|20d6r showed 
that A'p values for species like CO and CH in a single 
versus double sided calculation vary by up to a factor of 
two for an Ay between 0.2-5 mag. 

Our choice of explored ranges in physical parameters 
such as density, radiation field intensity, cosmic ray ion- 
ization rate, and stopping criterion (the physical thick- 
ness of our slab model) are determined by the need to 
compare our results with observations and with results 
from previous studies, and by typical diffuse cloud condi- 
tions. Diffuse clo uds generally have n n ranging from 10 
to 5000 cm~3 (Sn ow fc McCal]||2006f) . and therefore we 
vary nn from 10 to 1000 cm~'^, in increments of 1 dex. 
For simplicity, all calculations are perfor med at constant 
(depth-independent) density. We use the iDraing ()1978f ) 
radiation field in our calculations , which is also used in 
Ivan Dishoeck fc Blackl (|1988f) and lLe Petit et all (|2006D . 
We vary the far-UV intensity from /uv = 0.1 to 10 times 
the average value of the interstella r radiation fie ld, which 
equals 1.6 x 10"^ erg cm"^ s"^ (|Habing| [l968[ ) . also in 
increments of 1 dex. We also include the effect of cos- 
mic rays, for which we use a cosmic ray ionization rate 
C = 3 X 10~^^ cn i ^ s~^ . Hig her values o f C wer e found 
by iMcCall et all (l2003l) and iShaw et all (l2006h study- 
ingH3__toward ( Per and HD185418, respectively; by 
iLisztl ( 20031) . who inferred a higher ( based on analy- 
sis of_JiD_and_J|^ along a sample of sight lines; and 
bv lFederman et sH ()1996br ) based on cosmic ray-induced 
chemistry of OH toward o Per. Nevertheless, since our 
goal is to model global trends, we decided to use a value 
of C consistent with the average value o f C = 2.5 x 10""'^'^ 
cm^ s~^ determined by William s et al.l ((1998). We stop 
all calculations once iV(H2) = 2 x 10^^ cm"^, a value high 
enough to include all the diffuse cloud observational data 
in our sample. This stopping criterion corresponds to 
Av of 1-5 mag. We integrate molecular TV for all A^(H2) 
values up to the stopping criterion, t hus including all 
the phases of ISM clouds as classified by lSnow fc McCaOl 
(gOM). 

The thermal and ionization balance are both com- 
puted self-consistently. The temperature is computed 
from energy conse rvation consisting of a host of micro- 
physical processes (jFerland et al.lll998l : lAbel et al.ll2005l : 
iRolhg et al.ll2007t ). All atomic photo-processes are cal- 
culated by integrating the product of the incident radi- 
ation field intensity over the cross section for the photo- 
interaction rate. We also integrate the cross section for 
photodissociation of H2 , using a d etailed H2 model incor- 
porated into Cloudy (Shaw et al . 2005). For CO, we use 
the sh ielding function described in van Di shoeck fc BlacM 
(I1988D and iHollenbach. Takahashi. fc Tielend (|1991h . 
ILe Petit et al.l (|2006D show the predicted iVfCO) is about 
a factor of two smaller when using a shielding function 
versus an exact treatment of CO photodissociation in 
diffuse clouds. This is small, however, when compared 
to the increases by a factor of 100 that nonequilibrium 



CH+ formation can c ontribute to the formation of CO in 
diffuse environments (jZsargo fc Federmanll2003l and this 
work) . 

Our assumed gas and dust abundances are consistent 
with average ISM values. For the gas phase, we in- 
clude the 30 lightest elements. The abundance rela- 
tive to hydrogen for e ach species is an average of the 
abundance taken fr om 'Co wie fc Songailal (|1986l ) and C 
Oph (Sav age fc Semb ach 1991). The only exceptions are 
C/H and 0/H. For C/H w e use the value determined by 
ISavage fc Sembach! (|1996D without averaging, while fo r 
0/H we use the value determined by lMever et al.l (|1998f l. 
Some of the more important abundances by number are: 
He/H = 0.098, C/H = 1.3 x 10"'*, 0/H = 3.2 x 10-^ 
N/H = 8 X 10-5, Ne/H = 1.2 x lO""*, Si/H = 3.2 x 
10-^ S/H = 3.2 X 10-^ and Cl/H = 1 x 10"^ 

Our network includes all known important chemical 
channels leading to CO and CH"*" formation (Fig. 20). 
The chemical reaction network consists of approximately 
1200 reactions involving 89 molecules made up of H, He, 
C, N, O, Si, S, and CI. A complete list of molecules and 
reactions, along with rates, can be found on the Cloudy 
website^ . Most of the rate coefficients come from the 
UMI ST database (ILe Teuff et al.l l2000t iWoodall et all 
|2007( ). although there are a few exceptions. For the im- 
portant C+(OH,H)CO+ reaction, we use a tem perature- 
depeii dent rate based on the data of iDube rnct et aD 
()1992| ). with an equation derived in lAbel et al., (,20051). 
H2 is known to form primarily through catalysis on grain 
surfaces, and we compute the rate of H2 formation using 
the temperature- and m aterial-dependent rates given in 
ICazaux fc TielensI (|2002D . 

The most important aspect to our calculations is the 
modeling of nonequilibrium chemistry in order to sim- 
ulate the formation of CH+ and its trickle-down ef- 
fects_on CO. To t his eii d, we use the method given 
in Fe derman et al.l (|1996a[ l that incorporates a coupling 
between the ions and neutrals. The physical model 
for the nonequilibrium chemistry involves Alfven waves 
that, upon entering the cloud, dissipate over some phys- 
ical scale, as described in § 6.2. We model this ef- 
fect by reducing the coupling by one-third for N(H2) > 
4 X 10^° cm"^. This roughly corresponds to the tran- 
sition from the Diffuse At omic to Diffuse Molecular 
phase (Snow fc McCal]||2006D . and effectively "turns off" 
nonequilibrimn effects for N{H2) > 4 x 10^" cm"^. 
The coupling in terms of Tcs depends most critically 
on Auturb, which is the turbulent velocity of the gas. 
Therefore, we study the effects of At;turb on model pre- 
dictions by using three different values for it, 2.0, 3.3, 
and 4.0 km s"^. These values for Auturb are physically 
motivated from a number of con siderations. A t ypical 
CH+ linewidth is ~2.5 km s"^ (ICrawford et al.lll994 : 



iCrane et al.lll995tlCrawfo"rdlll995l : IPaIi et al.ll2004f ). The 
character of Awturb in the nonequilibrium chemistry is 
3-D; therefore if the CH+ linewidth is completely de- 
scribed by a 1-D turbulence driven by Alfven waves, 
Aftu rb = V^x linewidth ~ 4.3 km s~^ ([Heiles fc Trolandl 
l2005f ). These authors found that a typical 1-D value for 
Awturb ^ 1.2-1.3 km s"^ in the CNM, corresponding to 
2.1-2.3 km s~^ for 3-D turbulence. The average of the 
CH+ and CNM turbulence is ~3.3 km s^^. For all com- 
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binations of nn and /uv considered, we used this average 
value for Awturb (3.3 km s^^), and only compute models 
with Awturb = 2.0 or 4.0 km s^^ for nn — 100 cm~^ and 
-^uv — 1. More details about the nonequihbrium CH+ 
chemistry follow. 

6.2. Forming and Modeling CH^ 

Despite being one of the earliest molecules detected 
in the ISM, the formation of CH+ in the diffuse ISM 
remains one of the biggest challenges in astrochemistry. 
The fundamental issue is that equilibrium chemical mod- 
els under-predict iV(CH+) by 3-4 orders of magnitude. 
The main problem in the formation of CH"*" is the pri- 
mary formation channel leading to CH+ in the reaction 

C+(H2,H)CH+. (1) 

This reaction is highly endothe rmic, with a rate of 1 x 
10-i°exp(-4640/rkin) cm^ s'^ (jFederman et al.lll996aD . 
One way around the endothermicity of (1) is to have 
0+ react with excited II2 (II2*), which reduces or elim- 
inates the exponential temperature dependence in the 
rate constant. As was me ntioned in § 3.4, observations by 
iLambert fc Danki ()1986[ ) provided correlations between 
A^(CH+) and 7V(H2*) for J = 3 and 5, but our sample of 
sight lines does not show any similar correlations, a fact 
that could be the result of our much narrower range of 
examined A^(Il2*) values. In any case, models of CH+ 
chemistry, which include the formation process via II2*, 
still do not reproduc e the observed CH+ in diffuse clouds 
(|Garrod et al.l[2003l ) . CII+ can also form through the ra- 
diative association of 0+ and H. The rate for this reac- 
tion is 1.7 X 10~^^ cm'^ s~^, which exceeds the rate for 
reaction (1) for temperatures lower than 300 K. CII+ is 
easily destroyed, however, through reactions with H and 
H2. Therefore, the only way to efficiently produce CII+ 
to observed levels is to increase the temperature above 
the value used in models of equilibrium chemistry. 

It is generally agreed that nonequihbrium chemistry is 
the key to solving the CII+ abundance problem in dif- 
fuse clouds. However, the exact physical mechanism pro- 
ducing CII+ is still unclear. H ydrodynamic or Magne- 
tohydrodynamic sho ck models (jElitzur fc Watsonlll978l : 
iDraine fc Katd[T986l ) generate large amounts of CH+ by 
heating the gas to where CH+ efficiently forms by equa- 
tion (1). However, the lack of velocity differen ces be- 
tween CH and CH+ (Gre del et a l. 1993; Fedcrma n et al.l 
Il996al Il997b( ) argues agai nst shocks, w hile excitation 
analysis of interstellar C2 ()Grede]iri999f ) suggests CH+ 
production occurs in r egions where the gas temperature 
is 50-100 K. Recentlv.'Lesaf fre et all ()2007[ ) modeled the 
effects of turbulent diffusion on diffuse cloud chemistry, 
determining that this mechanism can increase the CH+ 
abundance by up to an order of magnitude, which is still 
~2 orders of magnitude lower than observed. 

It has been suggested that CH+-formation is driven by 
non-Maxwellian ve locity distributions of H2 and/or C+ 
([Gredel et al.lll993r ). One possible solution to the prob- 
lem of forming large quantitie s of CH+ in cold (7kin_< 
100 K) regions is discussed in iFederman et all ()1996aD . 
In this work, the authors propose Alfven waves entering 
a diffuse cloud from the intercloud medium are coupled 
to the cold gas through the Lorentz force (for ions) and 
collisions with the ions (for neutral atoms/molecules). 
The coupling results in significant nonthermal motion of 



the gas along the physical extent over which the MHD 
waves do not dissipate, consisting of a boundary layer 
on the cloud-intercloud surface. As a result, an effective 
temperature can be defined that characterizes the reac- 
tion between t wo specie s undergoing nonthermal motions 
(jFlower et al.llT985 : Federman et al.lll996a) 



Tkin + -^(Awturb)^. 
6K 



(2) 



In this equation, k is the Boltzmann constant, ^ is the 
reduced mass of the system, and Auturb is assumed to 
equal the Alfven speed. For turbulent velocities consis- 
tent with the observed linewidths of CH+, Toff is large 
enough to significantly increase the reaction rate of equa- 
tion (1), increasing 7V(CH+) to values consistent with 
observation. 

While this physical process is not the only possible ex- 
planation for the observed CH+ abundance, this method 
does have several important characteristics. One is that 
it allows for the formation of CH+ without heating the 
gas to temperatures inconsistent with the observed level 
of molecular excitation. This mechanism also explains 
the lack of OH toward f Per as a result of ion-neutral 
decoupling, where the magnetic field was coupled to the 
ions but not the neutrals. Thus Teff increased the re- 
action rate of ion-neutral reactions such as equation (1), 
but not neutral-neutral reactions such as 0(H2,H)0H, an 
important pathway to OH production at high tempera- 
tures. Finally, Alfven wave propagation and dissipation 
in a cold diffuse cloud is a relatively simple way to model 
nonequihbrium effects in a calculation designed to model 
equilibrium chemistry. All one needs to do is compute 
Tcff for each reaction using equation (2) and replace T 
with Toff when calculating the rate coefficient. 

6.3. Effects of CH+ on Other Molecules 

The regions where CH+ forms also contain signif- 
icant quantities of other molecules. This conclusion 
is independent of the ac tual physical p r ocesses con- 
troUing CH"*" formation. IFederman et al.l (|1997b[ ) and 
IZsargo fc FedermanI (j2003l ) estimated the contribution to 
the formation of CH and CO due to equilibrium processes 
alone (i.e. due to regions that do not form CH+) using a 
simple chemical model of a diffuse cloud. These studies 
found that most CH and CO (over 90% in many cases) 
could not be explained through equilibrium processes. 
The conclusion is that CH and CO in low-density (nn < 
100 cni^"^) sight lines form in regions where nonequihb- 
rium processes dominate the chemistry. So the same 
physical process that controls CH"*" formation is also 
likely to contribute to the formation of these molecules. 

Almost all explanations to account for CH+ involve in- 
creasing the temperature in order to activate the forma- 
tion channel given by equation (1). However, increasing 
the temperature also increases the rates of other reac- 
tions, leading to increased formation of certain molecules. 
One example, OH, has already been mentioned. Form- 
ing CH+ via equation (1) also leads to increased for- 
mation of CH through the chain involving CHj and 
CHj^, as was given in § 5.2. Forming CH+ also leads to 
the formation of C0+ via the reaction CH+(0,H)CO+, 
which is then followed by these two CO-forming chan- 
nels: CO+(H,H+)CO and C0+(H2,H)HC0+ together 
with HCO+(e,H)CO. Finally, CO and CH are coupled 
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through the neutral-neutral reaction CH(0,H)CO. This 
last reaction is an efficient formation route of CO and 
destruction route for CH at the high effective temper- 
atures required for CH"'' formation. At higher density 
(riH > 100 cm-3), C+(OH,H)CO+, followed by either 
CO^-to-CO channel above, becomes the primary route 
for CO formation. Regardless of nn, photodissociation 

is the primary destruction process fo r CO. 

Calculatio ns mad e by [ran Dishoe ck fc Blaclq (|1988f) . 
IWarin et all (1996"), and Le Petit et al." (2006) did 
not consider n oncquilibrium effe cts. The mo d els o f 
Ivan Dishoeckfc Black ( 1988l}_aiid [Le Petit et all (j2006l ) 
were used bv ISonnentrucker et al.l (|2007l ) to show that, 
if TT-H is sufffciently high, the correlation between CO 
and H2 observed through UV absorption from Coperni- 
cus, lUE, FUSE, and HST can be reproduced (see also 
Fig. 7). However, since neither model considered CH+ 
formation nor noncquilibrium chemistry (although the 
Meudon PDR group did in the p ast take CH+ into ac - 
count through a shock model, see iLe Petit et all l2004f ) . 
the CO relationship with H2 is likely to be much differ- 
ent in a model that also reproduces trends in CH+ versus 
H2 and CH+ versus CO. Such modeling and comparisons 
with observed trends in H2, CO, and CH+ abundances 
are the goals of our analysis. 

6.4. Comparing Model Results to Observation 

The results of our calculations are shown in Figs. 21, 
22, and 23. Each one shows plots of Af(CH+) ver- 
sus A^(H2), A^(CH+) versus A^(CO), and A^(CO) versus 
A^(H2). Figure 21 shows the results for /yv — 1, log nn 
= 1, 2, and 3 (Auturb — 3.3 km s~^); Fig. 22 the results 
for log riH = 2, /yv = 0.1, 1, and 10 (Awturb = 3.3 km 
s^^); and Fig. 23 the results for log rin — 2, /yv = 1, 
and Awturb = 2, 3.3, and 4 km s~^. 

6.4.1. CH+ versus H2 

The general observed trend in this plot is that A^(CH+) 
appears to saturate around 2 xlO^'^ cm~^, after which 
increasing A^(H2) does not result in increased A^(CH+) 
(recall Fig. 10). This trend is likely due to CH+ forma- 
tion by a mechanism that is acting over only a portion 
of the cloud. Our model mimics this effect by reducing 
Teff for A^(H2) > 4 X 10^" cm^^. Once iV(H2) becomes 
greater than this limit, the combination of smaller Tcff 
(due to "turning off" the noncquilibrium chemistry) and 
increased destruction of CH+ through reactions with H2 
leads to a decreased CH+ density and hence a saturated 
A^(CH+). Figures. 21 and 22 show that for any value 
of A^(H2) the value of A^(CH"'") is inversely related to 
nn/Iw- As nn/Zuv increases, the depth at which hy- 
drogen becomes predomi nately molecular [/ ( H2) > 0.8] 
decreases (see Fig. 15 of iLe Petit et al.lf2006t ). while an 
increased H2 abundance increases the destruction rate of 
CH"*" through the formation of CH^ and, eventually, CH 
(§ 6.3). This effect is greatest for nn/Zuv = 1000 cm"'^ 
(corresponding to log nn = 3, /uv = 1 in Fig. 21 or log 
riR = 2, /uv = 0.1 in Fig. 22) because /(H2) is then 
nearly 1, wh ereas for smaller nH// TTV values /(H2) never 
exceeds 0.8 ()Le Petit et al.l [20061 ) . Our models predict 
(in Fig. 23) that A^(CH"'") increases as Auturb increases. 
This trend is easy to understand as a temperature effect. 
Increasing Awturb increases T^s, which increases the rate 
of reaction (1). 



We find that reasonable assumptions of model parame- 
ters can explain the observed distribution in the A^(CH"'") 
versus A^(H2) plot. If we limit ourselves to only the case 
where /uv = 1 (Fig. 21), then over half of the data 
points lie in the region between the log nn = 1 and 3 
lines. About 75% of the rest of the data fall above the 
log riH = 1 line, i.e., where tt-h < 10 cm~'^. When the 
effects of turbulence (Fig. 23) are considered, then essen- 
tially all the observations are consistent with a suitable 
combination of rifj and Auturb- 

6.4.2. CH+ versus CO 

A plot of log 7V(CH+) versus log A^(CO) shows an 
important observational trend that can be understood 
through our calculations. Initially, A^(CH"'") increases 
with A^(CO), as was found in Fig. 10. Once iV(CO) 
reaches lO^'' cm~^, A^(CH"'") no longer increases, but 
levels off at A^(CH+) ~ 2 x lO^^cm"^. For even larger 
A^(CO), A^(CH+) appears to decrease. This trend can be 
understood as reflecting variations in njj (or nn/Iuy, if 
/uv differs significantly from unity). As njj increases, the 
amount of CH+ decreases, while the amount of CO in- 
creases (see Fig. 21). Therefore, the observations are well 
characterized by variations in nn, since for log iV(CO) < 
14, where CH+ and CO are coupled through the CH+ + 
O reaction, nn = 10-100 cm~^, while regions of higher 
CO and lower CH+ have nn > 100 cm~^. This conclu- 
sion does depend somewhat on Awturb, since for Auturb 
= 2 km s-i, the amount of Af(CH+) per 7V(C0) falls 
off significantly. However, for regions with values of /uv 
consistent with the average interstellar far-UV radiation 
field, the ratio of iV(CH+)/A^(CO) combined with our 
models is a good diagnostic of the density and hence the 
importance of noncquilibrium effects. 

6.4.3. CO versus H2 

Our calculations for the variation in A^(CO) versus 
A^(H2) show several important re sults. Comparmg the 
results of the equilibrium models of lLe Petit et al.l ()2006D 
with the log nn = 2; /uv = 1 calculations, we find 
that including CH+ chemistry can increase A^(CO) by 
a factor of 50-100 at low column densities. Such a dra- 
matic increase in CO cannot be attributed to a more 
rigorous treatment of the CO dissociation rate or geom- 
etry effects , both of which wo uld enter at the factor- 
of-2 level ()Le Petit et al.l I2006D . Instead, this points 
to the chemistry of CH+ as essential for understand- 
ing the abundance of CO in diffuse clouds, especially 
for 7V(C0) < W^-W^ cm-2. Equilibrium models need 
high riH values between 100-10 00 cm~^ to match the 
obser ved A^(CO) versus iy(H2) ()van Dish oeck fc Blackl 
,1988: iSonnentrucker et al.ll2007| ). Our calculations show 
that the value of nn can be an order of magnitude lower. 
This supports the conclusion of IZsargo fc FedermanI 
(|2003[ ) that noncquilibrium chemistry is important to the 
formation of CO in diffuse environments. For larger den- 
sities (or large nYi_/ I\jy), CH+ chemistry is less important 
due to destruction of CH+ through the formation of CHJ 
(§ 6.3). For regions where A^(H2) exceeds the cutoff for 
noncquilibrium effects, modeled A^(CO) values remain 
constant as a result of the decline in Tcff and thus in CO 
formation rates. However, observed A^(CO) can readily 
reach values higher than 10^^ cm^^ for higher A^(H2). 
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This is because, for large nn/^Vj the CO photodisso- 
ciation rate is less effective in destroying CO while the 
reaction C+(OH,H)CO+ still forms CO efficiently. The 
effect of Toff on iV(CO) is easily seen in Fig. 23, where 
increasing Awturb by a factor of 2 leads to A^(CO) values 
higher by ~4 dex. Overall, about 35% of the observations 
require densities ranging from 100 to 1000 cni"'^, with 
the rest requiring lower densities and hence the effect of 
nonequilibrium CH+ chemistry on CO production. 

From Fig. 21 we also see that modeled nn values do 
increase both along and across the (variable) slope of the 
CO versus H2 distribution. This confirms the results in § 
4.3, where these trends were derived qualitatively based 
on observed CN/CH+ ratios. The range of na between 
10 and 1000 cm^'^ from Cloudy is in good agreement with 
the analytic results from the CN chemistry in § 5. For 
sight lines without detected CN, analytic CH+ chemistry 
yielded very low values for nn , of which some 30% were 
between 3 and 30 cm~^ and thus in agreement with the 
lowest numerical values in Fig. 21. It is evident from 
the Cloudy results that CH+ resides in regions of lower 
gas density, whereas higher values of Hh correlate well 
with increased A^(CO). This is in excellent agreement 
with our pointing out in §§ 4 and 5 that higher-density 
gas is associated both with CN and CO, as well as with 
our empirical finding that uh is revealed by the observed 
ratio iV(CN)/iV(CH+). 

7. DISCUSSION 
7.1. The Ratio of CO to H2 

In this study we showed (Fig. 6) that the power-law 
relationship, N(CO) ex [N{H2)]^ , is observed to behave 
differently in two density regimes that control the pro- 
duction route for CO. For sight lines below the break at 
log 7V(H2, CO) = (20.4, 14.1), the relationship has B = 
1.5 ± 0.2, while above the break for higher-TV, higher-nn 
sight lines, it becomes steeper with B = 3.1 ± 0.7. The 
higher value for B is consistent with CO photochemical 
predictions (jvan Dishoeck fc Blackl [l988f ) of the transi- 
tion region between the diffuse and dark cloud regimes 
(panel b of Fig. 6 and Fig. 7), where UV shielding 
plays an important role (9 « 0.1) for iV(CO) « 10^^ 
and A^(Il2) w 10^^ cm^^. Throughout the plot the ver- 
tical dispersion in log iV(CO) has a full width of ±1.0, 
a range that is much larger than what is expected from 
observational uncertainties alone. This intrinsic disper- 
sion is influenced by the value of /uv/'^h, as can be seen 
in Figs. 7, 21, and 22. For /uv that is not far from 1, 
the width of the dispersion is reflecting the variability of 
nn, since we showed in § 4.3 (and Fig. 17) the changes in 
CN/CH+ both along and across the relationship of CO 
versus H2 (for a constant /uv = 1). A quantitative nu- 
merical confirmation from Cloudy was provided in Fig. 
21, where for log A^(H2) < 20.5 the dispersion in CO is 
seen to correspond to a range in rin between 10 and 100 
cm~^. 

The CO versus H2 relationship can be recast into 



N{CO) 

N{H2) 



[N{H2)] 



(B-l) 



which means that in the regime of diffuse clouds, the 
abundance of CO relative to H2 is not a constant fac- 
tor but varies between dependence on the square root 



of 7V(H2) (S - 1 w 0.5) for the lower values of N, 
and dependence on the second power (B — 1 « 2.1) of 
iV(H2) for clouds with log N > 20.4. This steeper depen- 
dence, though, gets shallower again once the transition 
into the dark cloud regime has occurred, just before CO 
uses up all the C atoms in the gas. In fact, for the as- 
sumed full conversion of C atoms into CO, the constant 
CO/H2 ratio means that the relative abundance of the 
two molecules is independent of Af(H2) for the highest 
N values. Such global variations have a bearing on Xqq 
whenever measurements include diffuse and translucent 
sight lines, since low values of CO/H2 for low N val- 
ues directly translate into higher values for Xqq. There 
is no doubt that this X-factor is dependent on physi- 
cal conditions that affect the abundance of CO in diffuse 
molecular clouds. However, since our sample involves es- 
sentially local clouds, it is not relevant to the issue of 
variations in Xq q over Galactic sca les, where metallicity 
can play a role (jStrong et al.ll2003 ). 

Our range of log CO/H2 ratios shows values between 
-7.58 and -4.68, with the single exception of HD 200775 
having a value of —3.88. The latter value is 47% of the 
value obtained from a full conversion of all carbon atoms 
into CO molecules, log (2 x C/H) = —3.55. As remarked 
earlier, full conversion is expected inside dark clouds, 
a regime associated with the F DR illuminated by HD 
200775. iFederman et~al] (|l98tf ) presented log CO/H2 
values between —7.37 and —5.30, i.e., overall lower val- 
ues that were the result of their sample of sight lines 
with lower iV( CO). For th e small samples toward Cep 
0B2 and 0B3 iPan et all ([2005) obtained values from 
—6.31 to —4.85 and from —6.42 to —5.95, respectively, 
with the f ormer range c learly incl uding higher-iV(CO) 
sight l ines. iBurgh et al.l ()2007t ) and ISonnentrucker et "aLl 
(|2007f ) presented restricted samples that ranged from 
—7.00 to —4.74 and from —6.56 to —4.56, respectively. 
The results from the smaller samples show good agree- 
ment with ours, albeit their ranges are narrower, as ex- 
pected. 

7.2. Connections to Molecular Clouds 

We showed the correspondence between significant 
molecular absorption and the presence of molecular 
clouds see n in emission for a number of directions in 
the past. iGredel et al.l (JT992) mapped the high lati- 
tude cloud respon sible for the absorption seen toward 
HD 210121, while 'Gredel et al.' (l9 9j) mapped t he CO 
emission around HD 154368. iFederman et al.l (|1994l ) 
indicated the sight lines probing molecular clouds as- 
sociated with sta rs in Ta urus, Ophiuchus, and Cep 
0B3, while Wann ier et al.l ([l999) did the same for the 
dark cloud B5 a nd stars in Per 0B2. Most recently, 
iPan et all (|2005D examined the correspondence between 
CO cloudlets seen in emission and stars in Cep 0B2. 

We can do the same for additional sight lines from 
the current survey. The most extensive sets of measure- 
ments probe molecular clouds in Chamaclcon (HD 93237, 
HD 94454, HD 96675, HD 99872, and HD 102065), the 
Southern Coalsack (HD 106943, HD 108002, HD 108639, 
HD 110434, HD 114886, HD 115071, and HD 115455), 
and Lupus (HD 137595, HD 140037, HD 144965, 
HP 147683). From the emission maps compiled by 
lAndersson et al.l (|2002i) . an interesting trend is discerned. 
Only directions with significant iV(CO) values (greater 
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than 10^'' cm"2) and CN (- 10^^ cm"^ when avail- 
able) lie within the CO contours. These are HD 96675, 
which probes the Cham I cloud, and HD 144965 and 
HD 147683, which pass through a cloud in Lupus. A par- 
ticularly interesting sight line for future study is toward 
HD 147683, where iV(CO) is about lO^^ cni~2_ Unfortu- 
nately, no CN data exist at the present time. 

With the aid of the SIMBAD site at the Centre 
de Donnees astronomiques de Stasbourg, we found 
other likely associations based on the similarity in 
I'hSR- The direction toward HD 30122 appears to 
be p robing the envelope of L1538 seen in CO emis- 
sion (lUngerechts fc ThaddeusI Il987l ) . The gas toward 
HD 36841 in the Ori OBI association may be related 
to that seen in emission from the reflection nebula 
IC 423 (iMaddalena et aHfl986') in the dark cloud L913 
([Clemens fc Barvainig|ll988l ). For all other sight lines, no 
clear correspondence could be found. 

7.3. Further Chemical Considerations 

Having applied the analytical expressions to extract 
gas densities from chemical schemes involving CN and 
CH+ to numerous sight lines, we now have a clearer un- 
derstanding of their limitations. While higher densities 
are found for CN-rich directions, as also found from our 
more comprehensive models, the correspondence between 
density and points on plots oi N{CO) versus A^(H2), etc. 
is rather weak (Fig. 17). The relationships involving 
N{CN)/N{CR) are s tronger (see left panel of Fig. 16 
and the discussion in lPan et al.r 2005'). The best corre- 
spondence is seen when A^(C0)/A^(H2) is plotted against 
N{CN)/N{C}1+) in the right panel of Fig. 16. This 
arises because CN o nly probes denser diffuse gas (e.g., 
ICardein et all fl99l iPanet all 12005). in which CH+ is 
more likely to be destroyed by H2. 

The limitations involve a number of factors. For a 
iven velocity compo nent, the amount of CH+-like CH 
Lambert et al.|[r990r ) is not easily obtained; the disper- 
sion in the relationship between A^(CH) and iV(CH+) for 
directions withou t detectable amounts of CN is too large 
(jPan et al.ll2005( l. Most sight lines, however, reveal ab- 
sorption from all three molecules at a given velocity (e.g.. 
Table 3). As for the CH+ chemistry, the amount of ma- 
terial along a line of sight and the strength of the local 
interstellar radiation field are not well known, especially 
for stars greater than a kpc away. More comprehensive 
models, combining the synthesis of CH+ and CN, are 
needed for the next level of understanding. Our goal is 
to apply models based on Cloudy to this problem. 

In this work, we presented a series of Cloudy-based 
calculations of diffuse cloud conditions that simultane- 
ously reproduce the observed H2, CH+, and CO abun- 
dances in these environments. Diffuse sight lines with 
A''(CO) < 10^''-10^^ cm^^ are well characterized by re- 
gions with nn/IuY < 100 cm~^, but only if the effects 
of CH^ are taken into account. Without the effects of 
nonequilibrium CH+ chemistry, equilibrium calculations 
predict too little CO per H2. This result appears ro- 
bust to uncertainties in the C/H abundance or the CO 
photodissociation rate. Furthermore, iV(CH+) increases 
with increasing N{GO), until N{CO) reaches lO^-lO^^ 
cm-2. For iV(CO) > lO^^^ cm^^ iV(CH+) appears to 
decrease for increasing A^(CO). Our models show this 
is likely a density effect, with A^(CO) increasing, and 



A^(CH+) decreasing, with increasing n-a. Last, the ob- 
served trend of A^(CH+) flattening out at a few times 
10^'^ cm~^ can be explained if the nonequilibrium chem- 
istry acts only over a certain physical size, s uch as the 
Alfvcn wave propagation formalism in Fcderman et al.l 
(1996a) . The observed scatter in iV(H2) with A^(CH+) 
is best explained through a combination of density ef- 
fects and the importance of nonequilibrium processes, 
parametrized in this work by At^turb- 

7.4. The Synoptic View 

The sight lines from our study have properties compa- 
rable to those inferred from both (a) H I self- absorption 
(HISA) clouds with weak or no CO emission and (b) CO- 
poor H2 gas revealed via 7-rays and far-infrared (FIR) 
emission. The former category has been investigated in 
recent 21-cm i&dio absor ytion surveys of the Galactic 
plane IjGibson et al.ll2005l ). showing that cold atomic hy- 
drogen gas is not necessarily associated with detections 
of CO emission. Since these clouds can be small (< 0.6 
pc) with n > 100 cm~^ and Tgpin < 50 K, the physical 
conditions in them are very similar to the clouds studied 
here, or in other words, they correspond to the intermedi- 
ate c ategory of diffuse molecular clouds (jSnow fc McCall 
l2006f) . One may, therefore, assume that despite CO non- 
detections via radio emission, CO is likely present in 
these clouds, albeit with low CO/H2 values determined 
by small values of uh/Ivv- Such clouds with low N{CO) 
should, i n principle, be d e tected via UV absorption. For 
example, iKlaassen et al.l (|2005| ) provide A^(CO) < 6 x 
10^^ cm~^ for a small HISA feature, an upper limit that 
excludes only the top 6% of our diffuse sight lines. 

In fact, an analy sis of a Galactic plane survey by 
iKavars et all (|2005l ) has determined that 60% of HISA 
features are associated with CO emission, w ith n ^ few x 
lOOcm-3 and 6 K < Tgpin < 41 K. Although lXavars et al] 
(|2005D suggest that these are "missing link" clouds be- 
tween the atomic and dense molecular varieties, we point 
out that this region in parameter space is occupied by dif- 
fuse molecular clouds. With H/H2 < 0.01 and CO/H2 
< 10~^ (Klaasscn et al. 2005) these are probably clouds 
that include the types of carbon-bearing molecular pho- 
tochemistries that were explored here. The inferred Tgpin 
values are lower than the Toi(H2) kinetic temperatures 
along diffuse sight lines, resembling more To2(C2) val- 
ues that are assoc iated with denser diffuse gas and the 
presence of ^^CO (jSheffer et al.|[2007t) . It wiU be interest- 
ing to see if these colder clouds are relate d to sight lines 
with v ery low ^^CO/^'^CO as observed bv lLiszt fc LucasI 
(|1998f ) using millimeter-wave absorption observations. 

A second category of CO-p oor H2 ga s has 
been revealed in FIR st udies (iReach et al.l 119941: 
Meverdierks fc HeithausenI 119961 : iDouglas fc Tavloj 



2007t). The survev of lReach et al.l (|l994D found infrared 
excess emission from cirrus clouds attributed to cold H2 
gas and dust, whereas only half of the clouds showed 
detectable levels of CO emission. When detected, CO 
was found to be subthermally excited with inferred nn ~ 
200 cm-3 and i^CO/^CO ratios between ~10 and >90, 
all indicating overlap with the parameter space of diffuse 
molecular clouds. CO-poor H2 gas has been inferred als o 
from Galactic surveys of 7 rays (jGrenier et al.l [2005) , 
which trace the gas content in the ISM and show an 
"excess" of 7 emission not associated with CO emission. 
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Indeed, iGrenier et al.1 ()2005D indicate that the CO-less 
gas is found around dense molecular clouds (that are 
detected via CO emission) and along bridges between 
cloud cores and atomic gas, precisely the sites where one 
would find gas known as diffuse molecular clouds. As is 
the case with the HISA clouds, we believe that CO is 
still there, albeit at low levels of abundance relative to 
H2 that are potentially observable via UV absorption 
but are not seen via current methods that detect CO in 



emission. 



Federman fc Willsonl ()1982D showed that a connection 
exists between diffuse molecular gas and dark clouds, 
namely, that there is good agreement in CH abun- 
dance and radial velocity between radio emission for dark 
clouds and optical absorption along nearby sight lines, 
the latter probing the outer envelopes of dark clouds. 
This CH conn ection has been exploited by Magnani 
and c olleagues (|Magnani fc Oneliolll995t iMagnani et al.l 
I2003D in deriving CO/H2 ratios for translucent sight 
lines, b ased on the tight correlation between CH and H2 
(jFederm an 1982; our § 3.3). Magnani at al. (2005) ob- 
served radio emission from CH along the Galactic plane, 
finding similarities to CO emission line profiles and in- 
ferring that the molecular gas has Uh < 1000 cm~'^. The 
material resembles the denser gas in our sample of dif- 
fuse molecular clouds, consistent with the presence of 
CO emission (see § 7.2). This connection is also related 
to OH emission that has been detected from intermedi- 
ate regions around de nser molecular CO-emitting clouds 
(jWannier et al.l [l993l ) and the lack of a correlation be- 
tween A^(^'^ CO) and 7V(0H) in ano ther sample of molec- 
ular clouds (|Goldsmithfc Lill2005D . 

8. CONCLUSIONS 

Our study of diffuse molecular clouds employed a new 
and extensive sample of sight lines with UV observations 
of CO and H2 to explore in detail the power-law relation- 
ship between the two species. The slope of log iV(CO) oc 
B X log A^(H2) was shown to require two components, 
one with B = 1.5 ± 0.2 for log A^(H2) < 20.4, and another 
B = 3.1 ± 0.7 for VUV sight lines with higher iV(H2). 
The break in slope arises from a change in CO produc- 
tion, with CH+ -I- O important at low iV(CO) and C+ 
+ OH at iV(CO) > IQi* cm-2. The ratio CO/H2 has 
a dependence on A^(H2) that results in an increase by 
~3.5 orders of magnitude over the range of log A^(H2) w 
19.5-22.0. Causes for variation in Xqo include (a) the 
nB_/I\jY ratio, which affects production and destruction, 
including self shielding, and (b) the metallicity of gas. 

Together with the CO and H2 we also analyzed new 
data for the carbon-bearing diatomic molecules CH, 
CH+, and CN (as well as C2) that are accessible through 
ground-based spectroscopy. The linear relationship be- 
tween 7V(CH) and 7V(H2) was confirmed again, both di- 
rectly using these two molecules, as well as indirectly by 
showing that the CO versus CH relationship follows that 
for CO versus H2 . After determining fitted relationships 
of both CO and CH versus H2 we were able to employ 
fit parameters in the prediction of iVp(H2) for three sight 
lines without H2 data. Analyzing A^(CH+) versus H2 and 



CO resulted in two more confirmations of the power-law 
break displayed by CO versus either H2 or CH, show- 
ing that this break separates the regime of low-density 
photochemistry from that involving high density. As for 
A^(CN), all our regression fits returned slopes with B < 
1.8, somewhat shallower than earlier reports. Since es- 
sentially all CN detections are along high-density sight 
lines, the absence of a detected break in slope for the 
(smaller) CN sample is not surprising. 

Many of the sight lines here are helping us to explore 
molecular environments that are associated with low njj ■ 
As such, these lines of sight probe regions where nonequi- 
librium CH"*" chemistry is dominating the production of 
CO, as confirmed by modeling with Cloudy. For those 
sight lines with higher nn it was possible to include (equi- 
librium) chemistry of CH, C2, and CN to predict molec- 
ular abundances and gas density. Such predictions were 
found to have tighter correlations when A^(C2) is part 
of the input into the chemical model. For the entire 
range of densities we showed that the empirical ratio 
7V(CN)/7V(CH+) is better suited than iV(CN)/iV(CH) 
as an indicator of the average nn along diffuse sight lines. 

We also considered rotational (excitation) tempera- 
tures in our modeling of CO and H2, showing that 
Toj(CO) does not vary for J = 1-3. On the other hand, 
Toj(H2) increases with J, with indistinguishable slopes 
between log Tqj and log Tqi for J = 2-4. Further analysis 
of the excitation of both molecules should help constrain 
the conditions in diffuse molecular clouds. 

As related in § 7.4, it is our understanding that the 
regime of low-iV(CO), Iow-tih diffuse molecular clouds is 
also sampled by a variety of non-UV observational meth- 
ods, which nonetheless result in a significant number of 
CO nondetections. Thus the true nature of diffuse molec- 
ular clouds is best revealed by synoptic knowledge ex- 
tracted from studies spanning the electromagnetic spec- 
trum from radio and FIR, through visible and UV, to 
gamma ray observations. 
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TABLE 1 
Stellar Data for Sight Lines with New Detections of CO^ 



Star 


Name 


Sp. 


V 


I 


b 


^lsr'' 


E(B-V) 


Ohelio'^ 


Refs'l 








(mag) 


(dog) 


(deg) 


(km s-1) 


(mag) 


(pc) 




BD+48 3437 




Bl lab 


8.73 


93.56 


-2.06 


14.4 


0.35 


6500 


1,10 


BD+53 2820 




BO IV:n 


9.96 


101.24 


-1.69 


12.8 


0.29 


4100 


2,2 


CPD-69 1743 




B2 Vn 


9.46 


303.71 


-7.35 


-8.0 


0.30 


4700 


3,3 


CPD-59 2603 


V572 Car 


07 V 


8.75 


287.59 


-0.69 


-11.6 


0.46 


2600 


3,3 


HD 12323 




09 


8.92 


132.91 


-5.87 


3.5 


0.23 


3600 


2,2 


HD 13268 




08 Vnn 


8.18 


133.96 


-4.99 


3.3 


0.36 


2400 


2,2 


HD 13745 


V354 Per 


09.7 H 


7.90 


134.58 


-4.96 


3.1 


0.46 


1600 


1,10 


HD 14434 




06.5 


8.59 


135.08 


-3.82 


3.1 


0.48 


4100 


2,2 


HD 15137 




09.5 V 


7.86 


137.46 


-7.58 


1.8 


0.35 


2700 


1,10 


HD 23180 


o Per 


Bl III 


3.86 


160.36 


-17.74 


-6.7 


0.30 


430 


2,2 


HD 23478 




B3IV 


6.69 


160.76 


-17.42 


-6.7 


0.28 


240 


1,11 


HD 24190 




B2 V 


7.45 


160.39 


-15.18 


-6.4 


0.30 


550 


1,10 


HD 24398 


C Per 


Bl lab 


2.88 


162.29 


-16.69 


-7.1 


0.34 


300 


4,11 


HD 30122 


HR 1512 


B5HI 


6.34 


176.62 


-14.03 


-10.9 


0.40 


220 


5,10 


HD 34078 


AE Aur 


09.5 Ve 


6.00 


172.08 


-2.26 


-8.4 


0.53 


450: 


1,11 


HD 36841 




08 


8.58 


204.26 


-17.22 


-17.1 


0.35 


1200 


1,10 


HD 37367 


HR 1924 


B2 IV-V 


5.99 


179.04 


-1.03 


-10.1 


0.42 


240 


2,2 


HD 37903 




B1.5 V 


7.84 


206.85 


-16.54 


-17.6 


0.32 


790 


4,10 


HD 43818 


11 /LU Gem 


BO II 


6.92 


188.49 


+3.87 


-11.9 


0.52 


1600 


2,2 


HD 58510 




Bl lab 


6.80 


235.52 


-2.47 


-18.8 


0.32 


4500 


1,10 


HD 63005 




07 


9.13 


242.47 


-0.93 


-18.5 


0.32 


5200 


1,10 


HD 91983 




O9.5/B0 lb: 8.58 


285.88 


+0.05 


-11.9 


0.29 


7000 


2,2 


HD 93205 


V560 Car 


03 V 


7.76 


287.57 


-0.71 


-11.6 


0.38 


3200 


2,2 


HD 93222 




08 


8.11 


287.74 


-1.02 


-11.6 


0.36 


1700 


2,2 


HD 93237 


DRCha 


B4IVe 


5.97 


297.18 


-18.39 


-10.9 


0.09 


310 


6,11 


HD 93840 




B1.5 lab 


7.79 


282.14 


+11.10 


-11.1 


0.16 


5700 


1,10 


HD 94454 




Bsm 


6.70 


295.69 


-14.73 


-11.0 


0.18 


330 


6,11 


HD 96675 




B6 IV/V 


7.6 


296.62 


-14.57 


-10.7 


0.30 


160 


2,11 


HD 99872 


HR 4425 


B3 V 


6.11 


296.69 


-10.62 


-10.4 


0.36 


230 


2,11 


HD 102065 




B2 V 


6.61 


300.03 


-18.00 


-10.1 


0.17 


170 


7,10 


HD 106943 




B7IV 


7.51 


298.96 


+1.14 


-8.3 


0.15 


500 


1,10 


HD 108002 




B2 la/ab 


6.95 


300.16 


-2.48 


-8.4 


0.32 


3400 


1,10 


HD 108610 




B3 IV/V 


6.92 


300.28 


+0.88 


-7.9 


0.15 


380 


1,11 


HD 108639 




Bl HI 


7.81 


300.22 


+1.95 


-7.8 


0.37 


110 


1,10 


HD 110434 




B8/9 HI 


7.55 


302.07 


-3.60 


-8.0 


0.05 


370: 


6,11 


HD 112999 


V946 Cen 


B6 Illn 


7.38 


304.17 


+2.18 


-6.6 


0.23 


340 


1,11 


HD 114886 




09 V 


6.89 


305.52 


-0.83 


-6.6 


0.40 


1000 


8,10 


HD 115071 


V961 Cen 


09.5 V 


7.97 


305.76 


+0.15 


-6.4 


0.53 


1200 


1,10 


HD 115455 




07.5 HI 


7.97 


306.06 


+0.22 


-6.3 


0.49 


2000 


6,10 


HD 116852 




09IH 


8.49 


304.88 


-16.13 


-8.6 


0.21 


4800 


2,2 


HD 122879 


HR 5281 


BO la 


6.43 


312.26 


+1.79 


-4.2 


0.36 


2300 


2,2 


HD 124314 




07 


6.64 


312.67 


-0.42 


-4.4 


0.53 


1100 


3,3 


HD 137595 




B3 Vn 


7.50 


336.72 


+18.86 


5.7 


0.25 


400 


1,10 


HD 140037 




B5III 


7.48 


340.15 


+18.04 


6.6 


0.09 


270: 


6,11 


HD 144965 




B3 Vne 


7.12 


339.04 


+08.42 


5.2 


0.35 


290 


6,10 


HD 147683 


V760 Sco 


B4 V 


7.05 


344.86 


+10.09 


7.2 


0.39 


280 


1,10 


HD 147888 


pOph D 


B3/B 4V 


6.78 


353.65 


+17.71 


10.5 


0.51 


140 


2,11 


HD 152590 




07.5 V 


8.48 


344.84 


+1.83 


6.2 


0.48 


1800 


1,10 


HD 152723 




07/08 


7.31 


344.81 


+1.61 


6.1 


0.42 


1600 


2,2 


HD 157857 




07e 


7.81 


12.97 


+13.31 


14.9 


0.43 


1900 


2,2 


HD 163758 




06.5 


7.32 


355.36 


-6.10 


8.1 


0.33 


2600 


4,10 


HD 185418 




B0.5 V 


7.52 


53.60 


-2.17 


18.1 


0.50 


910 


2,2 


HD 190918 


V1676 Cyg 


WN 


6.81 


72.65 


+2.07 


18.0 


0.45 


2300 


1,10 


HD 192035 


RXCyg 


BO Illn 


8.22 


83.33 


+7.76 


17.3 


0.37 


2800 


1,10 


HD 192639 




08e 


7.11 


74.90 


+1.48 


17.7 


0.62 


1600 


4,10 


HD 195965 




BO V 


6.98 


85.71 


+5.00 


16.7 


0.25 


790 


1,10 


HD 198781 


HR 7993 


B0.5 V 


6.46 


99.94 


+12.61 


14.7 


0.35 


730 


2,2 


HD 200775 


V380 Cep 


B2 Ve 


7.42 


104.06 


+14.19 


13.9 


0.57 


430: 


9,11 


HD 203532 


HR 8176 


B3IV 


6.36 


309.46 


-31.74 


-8.6 


0.28 


250 


2,11 


HD 208905 




Bl Vp 


7.01 


103.53 


+5.17 


13.1 


0.37 


790 


1,10 


HD 209481 


14/LZ Cep 


09 V 


5.55 


102.01 


+2.18 


13.1 


0.37 


690 


1,10 


HD 209975 


19 Cep 


09Ib 


5.11 


104.87 


+5.39 


12.8 


0.34 


1300 


2,2 


HD 210121 




B9 


7.69 


56.88 


-44.46 


7.9 


0.31 


210 


2,11 


HD 210809 




09Ib 


7.56 


99.85 


-3.13 


13.0 


0.31 


4000 


2,2 


HD 220057 


NSV 14513 


B2IV 


6.95 


112.13 


+0.21 


10.4 


0.23 


560 


2,2 


HD 303308 




03 V 


8.21 


287.59 


-0.61 


-11.6 


0.45 


3600 


2,2 


HD 308813 




09.5 V 


9.32 


294.79 


-1.61 


-10.0 


0.31 


2400 


1,10 
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^ Information from the SIMBAD database is included. Correction from heliocentric velocity to the LSR frame. ^ Distance 
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reference, or from a >4 a Hipparcos parallax from reference 10 as listed by Simbad, unless a >3 (T parallax was used and 
flagged with ":". First reference is for E{B — V), the second is for .Dhciio- 
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TABLE 2 
UV Data Sets for New CO Sight Lines 



Star 


_ff5T/STIS 




FUSE 
Data set 


Star 




ffSr/STIS 




FUSE 




Data sot 


Grating 


Slit 


Data set 


Grating 


Slit 


Data set 


BD+48 3437 


o6359s 


E140M 


0.2X0.2 


P10184 


HD 110434 


o61j0b 


E140H 


0.1X0.03 


A12019 


BD+53 2820 


o6359q 


E140M 


0.2X0.2 


P12232 


HD 112999 


o61j0c 


E140H 


0.1X0.03 


A12020 


CPD-69 1743 


063566 


E140M 


0.2X0.2 


P10137 


HD 114886 


o61j0d 


E140H 


0.1X0.03 


A12018 


CPD-59 2603 


o40p01 


E140H 


0.2X0.09 


P12215 


HD 115071 


o61j0c 


E140H 


0.2X0.09 


G93215 




o4qx03 


E140H 


0.2X0.09 




HD 115455 


o61j0f 


E140H 


0.1X0.03 


A12007 


HD 12323 


O63505 


E140M 


0.2X0.2 


P10202 


HD 116852 


063571 


E140H 


0.2X0.2 


P10138 


HD 13268 


O63506 


E140M 


0.2X0.2 


P10203 


HD 122879 


o61z57 


E140H 


0.2X0.2 


B07105 


HD 13745 


o61z05 


E140M 


0.2X0.2 


P10204 


HD 124314 


O54307 


E140H 


0.1X0.03 


P10262 


HD 14434 


O63508 


E140M 


0.2X0.2 


P10205 




o61z58 


E140H 


0.2X0.2 




HD 15137 


o61z06 


E140H 


0.2X0.2 


P10206 


HD 137595 


o61j03 


E140H 


0.2X0.09 


A12012 


HD 23180 


O64801-4 


E140H 


0.2X0.05 




HD 140037 


o61j04 


E140H 


0.1X0.03 


A12015 


HD 23478 


06IJOI 


E140H 


0.1X0.03 


A12002 


HD 144965 


o61j05 


E140H 


0.1X0.03 


A12016 


HD 24190 


06IJO2 


E140H 


0.1X0.03 


A12001 


HD 147683 


06IJO6 


E140H 


0.2X0.09 


A12009 


HD 24398 


O64810-11 


E140H 


0.2X0.05 




HD 147888 


o59s05 


E140H 


0.2X0.09 


P11615 


HD 30122 


o5c065 


E140H 


0.2X0.2 


Q20103 


HD 152590 


o61z67 


E140M 


0.2X0.2 


B07106 


HD 36841 


063516 


E140M 


0.2X0.2 




HD 152723 


063586 


E140H 


0.2X0.2 


P10271 


HD 37367 


o5c013 


E140H 


0.2X0.2 


B07102 


HD 157857 


o5c04d 


E140H 


0.2X0.2 


P10275 


HD 37903 


o59s04 


E140H 


0.2X0.09 


P11606 


HD 163758 


063595 


E140H 


0.2X0.2 


P10159 


HD 43818 


o5c07i 


E140H 


0.2X0.2 




HD 185418 


o5c01q 


E140H 


0.2X0.2 




HD 58510 


O63530 


E140H 


0.2X0.2 


P10219 


HD 190918 


o6359j 


E140M 


0.2X0.2 


P10285 


HD 63005 


063531 


E140M 


0.2X0.2 


P10221 


HD 192035 


o6359k 


E140M 


0.2X0.2 


P10286 


HD 91983 


o5c08n 


E140H 


0.2X0.2 


B07104 


HD 192639 


o5c08t 


H140H 


0.2X0.2 




HD 93205 


o4qx01 


E140H 


0.2X0.09 


P10236 


HD 195965 


o6bg01 


E140H 


0.1X0.03 


P10288 


HD 93222 


o4qx02 


E140H 


0.2X0.09 


P10237 


HD 198781 


o5c049 


E140H 


0.2X0.2 


P23102 


HD 93237 


o61j0g 


E140H 


0.1X0.03 


A12010 


HD 200775 








A05101 


HD 93840 


063549 


E140H 


0.2X0.2 


P10127 


HD 203532 


o5cols 


E140H 


0.2X0.2 


B07108 


HD 94454 


o61j0h 


E140H 


0.1X0.03 


A12005 


HD 208905 








D01401 


HD 96675^ 


zl9w01=' 


GieOM'' 


0.25^ 


Q10102 


HD 209481 








D01402 


HD 99872 


o61j0i 


E140H 


0.1X0.03 


A12006 


HD 209975 








D01403 


HD 102065 


o4o001 


E140H 


0.2X0.09 


QlOlOl 


HD 210121 








P24901 


HD 106943 


o61j07 


E140H 


0.1X0.03 


A12011 


HD 210809 


o6359t 


E140M 


0.2X0.2 


P12231 


HD 108002 


06IJO8 


E140H 


0.1X0.03 


A12017 


HD 220057 


o5c01x 


E140H 


0.2X0.2 


Z90178 


HD 108610 


o61j09 


E140H 


0.1X0.03 


A12014 


HD 303308 


o4qx04 


E140H 


0.2X0.09 


P12216 


HD 108639 


o61j0a 


E140H 


0.2X0.09 


A12013 


HD 308813 


063559 


E140M 


0.2X0.2 


P12219 



The HST data set for HD 96675 is from the GHRS, not the STIS. 
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TABLE 3 

Newly Detected Cloud Components of Molecular Species 





i'lsr'' 


CO 




CN 




CH 




CH+ 




Star 


N" 


j,a 


N" 


^a 


N" 


^a 


ATb 


{,a 






(10") 




(1012) 




(1012) 




(1012) 




BD+4S3437 


-16.4 

-11.4 














4.35 
0.98 


1.8 
1.0 




-1.0 


0.23 


1.5 










2.15 


1.7 




5.2 






0.41 


0.5 


6.86 


3.2 


7.25 


1.6 


BD+532820 


-5.4 














2.15 


1.0 




0.9 


0.50 


0.5 






1.59 


0.5 








6.9 


0.27 


0.4 






2.49 


1.1 


4.26 


2.2 


CPD-691743 


0.8 


0.11 


1.4 










15.42 


1.4 


CPD-592603 


-5.7 


0.08 


1.5 






3.80 


1.3 








-2.8 


1.13 


0.7 






7.91 


0.5 


16.36 


2.9 


HD 12323 


-13.5 


2.53 


0.7 


0.63 


0.5 


1.97 


1.0 








-9.7 






0.73 


0.5 


2.37 


1.0 








-5.9 


0.89 


0.8 






1.98 


2.3 


5.52 


2.5 




0.3 














2.42 


1.5 


HD 13268 


-36.2 
-19.5 










0.94 


1.0 


3.1 


3.0 




-16.4 


0.16 


1.3 






1.13 


1.0 








-10.4 






0.43 


0.5 


2.06 


1.9 


3.5 


2.3 




-7.4 


1.26 


1.0 


0.94 


1.3 


5.29 


2.3 








-5.1 














4.8 


1.6 




-1.0 


0.14 


1.0 






2.15 


2.5 


3.5 


1.5 


HD 13745 


-43.9 


0.33 


0.7 


0.52 


1.1 


5.85 


2.3 


8.59 


2.0 




-18.1 


0.36 


0.5 






3.97 


1.9 


10.98 


2.8 




-10.3 


0.11 


1.3 










7.22 


3.0 




-4.3 


0.07 


0.3 
















0.0 














5.89 


3.0 


HD 14434 


-6.1 


0.43 


0.6 










8.52 


3.0 




-1.0 


1.83 


0.7 


0.39 


0.8 


9.24 


2.3 


14.94 


2.5 


HD 15137 


-13.4 


0.11 


0.7 






2.17 


2.5 


4.13 


2.9 




-7.4 










1.49 


1.0 


7.70 


3.0 




-4.0 


0.03 


0.6 
















-0.2 


0.19 


1.3 






2.50 


2.1 


2.26 


1.7 


HD 23180 


4.6 


0.81 


0.8 






7.01 


2.1 


1.58 


2.0 




7.3 


5.97 


1.0 


1.33 


1.7 


11.96 


1.5 


5.72 


2.0 


HD 23478 


4.1 


6.28 


1.0 


0.67 


1.1 


13.53 


1.8 


1.42 


1.2 




7.7 


1.77 


0.9 


0.43 


0.5 


4.72 


0.9 


1.40 


1.0 


HD 24190 


6.5 
9.7 


0.84 
0.06 


0.5 
1.0 


0.50 


1.8 










HD 24398 


6.8 


17.95 


0.8 


3.20 


1.0 


21.41 


1.6 


3.13 


2.3 


HD 30122 


4.2 














2.98 


2.9 




6.9 


7.04 


0.8 


0.87 


1.3 


15.66 


2.0 






HD 36841 


-4.6 
5.7 


0.46 
0.35 


0.5 
0.4 
















10.6 






0.78 


0.5 


9.98 


1.5 


5.74 


2.2 


HD 37367 


3.8 










2.21 


2.0 








6.2 


0.70 


1.5 






9.80 


2.0 


32.41 


2.2 


HD 43818 


-7.0 










3.22 


2.5 








-3.9 


0.74 


1.5 






3.23 


1.6 


3.61 


1.8 




1.2 










2.77 


1.6 


7.36 


3.0 




5.2 










2.03 


1.9 


4.14 


3.0 


HD 58510 


23.6 














1.77 


1.5 




26.7 










2.52 


2.0 


10.74 


2.1 




29.6 


0.18 


1.5 






2.53 


1.6 






HD 63005 


9.9 














2.67 


1.8 




14.3 


0.63 


0.9 


0.99 


1.1 


4.51 


2.2 


5.67 


2.5 




21.0 


0.39 


1.5 






4.46 


2.0 


5.69 


1.7 


HD 91983 


-14.0 


0.62 


0.5 










3.27 


1.2 


HD 93205 


-6.5 


0.06 


1.0 






1.33 


1.0 








-2.8 


0.08 


0.3 






1.31 


0.5 






HD 93222 


-6.3 


0.18 


1.5 






2.28 


2.1 






HD 93237== 


3.5 


0.25 


0.7 






1.2 


0.7 






HD 93840 


-7.0 


0.18 


0.8 






1.79 


0.5 


3.28 


2.7 


HD 94454<: 


3.6 


2.02 


1.3 






9.4 


1.6 






HD 96675 


4.1 


20.18 


0.9 


4.96 


0.5 


22.76 


1.1 


4.90 


2.0 


HD 99872<: 


3.2 


4.54 


0.7 






12.6 


1.7 


23.10 


1.9 


HD 102065 


1.0 










1.14 


0.5 


4.97 


1.9 




3.8 


0.49 


1.4 






6.03 


1.6 


5.78 


1.2 


HD 106943 


-0.8 


0.06 


0.8 














HD 108002== 


-10.1 


0.06 


0.4 
















1.6 


0.35 


0.8 






3.2 


1.5 






HD 108639 


-1.2 


0.15 


1.0 
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TABLE 3 

Newly Detected Cloud Components of Molecular Species (Cont.) 



CO 



Star 



''LSR 



(10") 



CN 

(1012) 



CH 



CH+ 



(1012) 



(1012) 



HD 110434 
HD 114886'! 



HD 115071^^ 
HD 115455= 

HD 116852 
HD 122879 

HD 124314 

HD 137595== 
HD 140037 
HD 144965== 
HD 147683== 



HD 148937 

HD 152590 
HD 152723 

HD 157857 

HD 163758 
HD 177989 

HD 185418= 
HD 190918 

HD 192035 

HD 192639*" 
HD 195965 

HD 198781 

HD 203532 

HD 210121 

HD 210809 

HD 220057 

HD 303308 

HD 308813 



-5.1 

-1.7 

-27.5 

-4.3 

-1.9 

-3.3 

-3.3 

-0.1 

0.6 

-26.1 

-2.5 

-23.6 

-2.1 

5.4 

1.1 

4.5 

5.0 

6.1 

11.8 

-13.5 

-10.4 

3.2 

-0.5 

5.3 

2.1 

7.7 

10.2 

-5.6 

0.0 

4.2 

-3.8 

2.1 

6.7 

10.1 

12.1 

6.8 

11.0 

2.1 

5.8 

18.3 

1.4 

5.6 

9.4 

7.3 

6.4 

10.2 

-0.1 

5.4 

5.3 

7.5 

-6.2 

-1.7 

-1.2 

2.8 

-1.8 

1.9 

-7.7 

-4.5 

-1.4 

-2.8 



0.04 
0.04 
0.09 
0.19 
0.13 
3.40 
1.18 



0.13 
0.24 
1.33 
0.77 
0.41 
19.35 
0.54 
88.55 
0.18 
1.07 
2.51 
0.23 
0.44 
0.15 
0.10 
0.05 
0.43 



0.24 

0.04 
0.45 
3.96 
4.53 
0.80 

0.05 
0.09 

10.76 
0.81 
0.60 
0.50 
0.70 
0.66 
15.87 
45.62 



0.10 
0.35 
0.07 
0.69 



0.6 
0.3 
0.6 
1.5 
1.5 
1.3 
1.0 



0.20 0.5 



1.5 
1.5 
0.8 
1.0 
1.4 
0.6 
0.3 
0.6 
1.0 
0.4 
0.7 
0.8 
0.5 
1.5 
1.5 
0.5 
1.5 



1.22 0.( 



0.9 

1.0 
0.8 
0.5 
0.7 
0.3 

0.3 
1.5 

1.0 
1.5 
1.0 
0.5 
1.0 
0.3 
0.4 
0.6 



0.23 1.5 
4.34 1.1 



1.5 
1.4 
0.3 
0.7 



13 







05 O.f 



13 0.5 



34 0.9 



.96 0.( 



4. 




6. 
13. 




5. 








8.2 


2.7 






17. 


3.6 


9.06 


1.4 






3.30 


1.0 


1.73 


0.5 


1.93 


1.0 






3.94 


2.2 


2.45 


1.8 


7.24 


2.3 


1.75 


0.9 


8.59 


2.4 


7.07 


1.7 


6.23 


2.6 


12.2 


3.5 


9.63 


2.5 


14.4 


1.6 


7.58 


1.9 


2.0 


<0.3 


11.08 


2.5 


15.2 


1.9 


5.59 


2.8 


5.0 


2.5 


2.53 


1.6 


7.59 


1.9 


8.45 


2.3 


3.91 


3.0 


10.69 


2.1 


4.06 


5.2 


2.65 


1.5 


4.20 


0.8 


3.63 


1.9 


3.70 


1.4 


4.78 


2.2 






3.37 


2.3 


5.02 


1.0 


7.88 


2.8 


2.83 


2.6 


8.27 


3.0 


2.25 


1.1 







8.6 
3.2 
1.74 

1.16 

2.84 

10.12 

3.53 

6.8 



13.17 
9.76 
14.91 
28.61 

1.61 
3.86 

8.22 
4.87 

5.40 

5.09 



1.1 
0.5 
1.5 

0.6 
1.0 
1.0 
1.0 
1.2 



1.5 
3.1 
1.0 
1.6 

2.0 
1.5 
1.3 
2.5 

2.9 

0.9 



1.42 



3.9 
8.6 
3.77 
5.51 
4.74 
3.61 
4.11 



3.34 
2.98 

5.05 
6.55 

7.51 
5.27 



1.0 



1.3 
2.1 
1.4 
2.0 
2.5 
2.0 
2.0 



1.7 
2.0 

2.8 
1.7 

1.7 
3.0 



Units arc km : 



-1 b 



Unit; 



s arc cm 



-2 c 



CH results arc from lAndersson et al1l2002t'! CH and CH+ 



jlts a re from[Gr cdcl igQU " CH and CH+ results arc from lSonnentrucker et al.ll2003 .^ CH 
are from lSonnentrucker et al1l2007t 
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TABLE 4 
New Total Molecular Column Densities" 









log N (cm 


-') 








Star 


H2 Ref 


CO Ref 


CH+ Ref 


CH 


Ref 


CN Ref 


BD+48 3437 


20.42 


13.36 


13.18 




12.79 




11.76 


BD+53 2820 


20.15 


13.89 


12.81 




12.61 




<11.97 


CPD-69 1743 


19.99 


13.08 


13.18 










CPD-59 2603 


20.15 


14.08 


13.20 




13.08 






HD 12323 


20.32 


14.53 


12.90 




12.63 




12.30 


HD 13268 


20.51 


14.20 


13.18 




13.04 




12.30 


HD 13745 


20.67 


13.94 


13.52 




12.99 




11.86 


HD 14434 


20.43 


14.36 


13.38 




12.96 




11.74 


HD 15137 


20.32 


13.52 


13.15 




12.79 




<11.62 


HD 23180 


20.61 1 


14.83 2 


12.84 




13.28 




12.18 


HD 23478 


20.57 


14.91 


12.32 




13.34 




12.04 


HD 24190 


20.38 


13.95 


13.18 




12.98 




11.88 


HD 24398 


20.67 


15.26 


12.45 




13.32 




12.51 


HD 30122 


20.70 


14.85 


12.48 




13.20 




12.08 


HD 34078 


20.88 


14.76 


13.84 




13.90 




12.52 


HD 36841 


b 


14.08 


12.76 




13.00 




12.04 


HD 37367 


20.61 


13.85 


13.51 




13.08 




<11.53 


HD 37903 


20.95 


13.69 


13.11 


3 


12.96 


3 


11.90 3 


HD 43818 


b 


13.87 


13.18 




13.04 




<11.79 


HD 58510 


20.23 


13.26 


13.08 




12.70 




<11.71 


HD 63005 


20.23 


14.00 


13.15 




12.95 




12.15 


HD 91983 


20.23 


13.79 


12.52 










HD 93205 


19.83 


13.15 


<12.20 




12.42 






HD 93222 


19.84 


13.26 






12.36 






HD 93237 


19.80 


13.40 






12.08 


4 




HD 93840 


19.28 


13.26 


12.52 




12.26 






HD 94454 


20.76 


14.30 






12.97 


4 




HD 96675 


20.86 


15.28 


12.69 




13.36 




12.77 


HD 99872 


20.52 


14.65 


13.36 




13.11 


4 


<11.79 


HD 102065 


20.56 


13.69 


13.04 




12.86 




<12.19 


HD 106943 


19.81 


12.76 






<12.49 


4 




HD 108002 


20.34 


13.66 






12.51 


4 




HD 108610 


19.86 








<12.56 


4 




HD 108639 


20.04 


13.18 






<12.38 


4 




HD 110434 


19.90 


12.94 






<12.23 


4 




HD 112999 


20.11 


<13.23 












HD 114886 


20.34 


13.61 


13.28 


5 


12.95 


5 




HD 115071 


20.69 


14.53 






12.91 


4 




HD 115455 


20.58 


14.08 


13.23 




13.23 


4 


<12.20 


HD 116852 


19.83 


13.30 


12.54 




12.23 






HD 122879 


20.36 


13.11 


13.08 




12.38 






HD 124314 


20.52 


14.20 


13.18 




12.94 






HD 137595 


20.62 


13.89 


13.26 




13.08 


4 


<11.90 


HD 140037 


19.34 


13.61 






<11.85 


4 




HD 144965 


20.79 


15.28 


12.88 




13.15 


4 




HD 147683 


20.74 


15.95 


13.28 




13.34 


4 




HD 147888 


20.58 


15.28 


12.88 


6 


13.34 


6 


12.32 6 


HD 152590 


20.51 


13.77 


13.28 




13.08 






HD 152723 


20.30 


13.76 


13.04 




13.08 






HD 157857 


20.69 


14.08 


13.30 




12.89 




<12.06 


HD 163758 


19.85 


13.38 


12.15 




12.34 






HD 190918 


19.95 


13.18 


13.15 




12.46 




<11.53 


HD 192035 


20.68 


15.15 


12.89 




13.20 




12.71 


HD 192639 


20.75 7 


13.78 


13.61 


7 


13.45 


7 


<11.85 8 


HD 195965 


20.34 


14.08 












HD 198781 


20.56 


15.23 


12.52 




13.11 




12.46 


HD 200775 


21.15 


17.29 


12.97 


9 


13.51 


9 


13.08 9 


HD 203532 


20.70 


15.66 


12.48 




13.40 






HD 208905 


20.43 


14.62 


12.78 


6 


12.73 


6 




HD 209481 


20.54 


14.60 


12.72 


6 


12.83 


6 




HD 209975 


20.15 


14.04 


13.38 


6 


12.93 


6 




HD 210121 


20.86 


15.83 10 


13.08 




13.46 




13.20 


HD 210809 


20.00 


13.36 


12.88 




12.74 




<11.81 


HD 220057 


20.34 


14.63 


12.87 




13.11 




12.11 


HD 303308 


20.15 


13.72 






12.73 






HD 308813 


20.30 


13.84 






12.71 
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(2) IShcffcr ot al.l 12007J; (3) 


'Knauth et al.l 120011: 
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iThorburn et al.ll2003l; f 9) iFederman et ahi,1997ai: ClO) .Sonnent 


2004; (7) Rachford ct al.l 120021; (8) 
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" Supplemented by literature values from listed references. HD 36841 and HD 43818 
predicted here to have log Af(H2) = 20.4 ± 0.1, see § 4.4. 
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TABLE 5 
Power-Law Fits of Column Density Correlations^ 



CL 



B 



og Af(CO) log Af(H2) 105 0.834 >99.99% -24.4 ± 3.1 1.89 ± 0.15 

og Af(CO) log N{'R2) < 20.4 50 0.734 >99.99% -15.8 ± 4.5 1.46 ± 0.23 

og Af(CO) log Af(H2) > 20.4 55 0.638 >99.99% -48.8 ± 15.0 3.07 ± 0.73 

og iV(CH) log Af(H2) 90 0.906 >99.99% -6.80 ± 1.50 0.97 ± 0.07 

og iV(CH) log Af(H2) < 20.4 36 0.799 >99.99% -5.87 ± 3.78 0.92 ± 0.19 

og Af(CH) log Af(H2) > 20.4 54 0.740 >99.99% -9.34 ± 3.90 1.09 ± 0.19 

og Af(CO) log Af(CH) 92 0.824 >99.99% -12.3 ± 2.7 2.05 ± 0.21 

og Af(CO) log Af(CH) < 13.0 42 0.624 >99.99% -5.30 ± 3.8 1.50 ± 0.30 

og 7V(CO) log 7V(CH) > 13.0 50 0.693 >99.99% -22.3 ± 11.3 2.80 ± 0.85 

log N(CU+) log Af(H2) 86 0.471 >99.99% 4.41 ± 2.06 0.42 ± 0.10 

log Af(CH+) log N(U2) < 20.3 26 0.637 99.95% -2.76 ± 4.30 0.78 ± 0.22 

IogAf(CH+) log Ar{H2) > 20.3 60 0.089 50% 10.1 ± 4.3 0.15 ± 0.21 

log Af(CH+) log N{CO) 88 0.120 73% 12.4 ± 0.7 0.04 ± 0.05 

log Af(CH+) log N(CO) < 14.1 41 0.648 >99.99% 6.83 ± 1.29 0.46 ± 0.10 

IogAf(CH+) log Ar{CO) > 14.1 47 0.268 93.1% 15.2 ± 1.0 -0.14 ± 0.07 

log Af(CN) log N{U2) 40 0.669 >99.99% -18.5 ± 7.5 1.49 ± 0.36 

log 7V(CN) log Af{H2) < 20.68 20 0.318 83% -9.0 ± 20.8 1.02 ± 1.02 

log Af(CN) log Af{H2) > 20.68 20 0.294 79% -11.5 ± 24.9 1.16 ± 1.19 

log Af{CO) log N{aSS) 42 0.836 >99.99% -2.81 ± 2.77 1.44 ± 0.23 

log Ar{CO) log Ar{CN) < 12.31 20 0.453 95.5% 2.69 ± 5.02 0.97 ± 0.42 

log Ar{CO) log Ar(CN) > 12.31 22 0.578 99.5% 3.38 ± 5.70 0.96 ± 0.45 



BCES(y|x) results from the lAkritas fc Bershadvlll996l code. 
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TABLE 6 
New Interstellar Excitation Temperatures for H2 and CO 



Star 



T„i(H2) 
(K) 



To2(H2) 

(K) 



703(H2) 

(K) 



To4(H2) 

(K) 



Toi(CO) To2(CO) To3(CO) 
(K) (K) (K) 



BD+48 3437 


83 


113. 


158. 


246 


BD+53 2820 


93 


120. 


176. 


244 


CPD-69 1743 


79 


102. 


143. 


213 


CPD-59 2603 


77 


95. 


142. 


217 


HD 12323 


82 


101. 


142. 


217 


HD 13268 


92 


120. 


167. 


245 


HD 13745 


66 


93. 


128. 


202 


HD 14434 


99 


129. 


166. 


247 


HD 15137 


104 


111. 


153. 


245 


HD 23478 


55 


79. 


101. 


171 


HD 24190 


66 


86. 


119. 


193 


HD 24398 










HD 24534 


54 


73. 


96. 


152 


HD 27778 


51 


78. 


103. 


152 


HD 30122 


61 


86. 


121. 


185 


HD 34078 


75 


92. 


128. 


206 


HD 36841 










HD 37367 


73 


82. 


112. 


185 


HD 37903 


64 


121. 


125. 


190 


HD 43818 










HD 58510 


90 


99. 


143. 


212 


HD 63005 


78 


91. 


129. 


188 


HD 91983 


61 


105. 


144. 


222 


HD 93205 


97 


118. 


167. 


241 


HD 93222 


69 


109. 


162. 


218 


HD 93237 


58 


85. 


111. 


135 


HD 93840 


54 


112. 


170. 


224 


HD 94454 


74 


83. 


106. 


167 


HD 96675 










HD 99872 


66 


94. 


114. 


179 


HD 102065 










HD 106943 


96 


108. 


142. 


214 


HD 108002 


77 


98. 


133. 


218 


HD 108610 


80 


106. 


138. 


208 


HD 108639 


88 


111. 


153. 


219 


HD 110434 


87 


105. 


144. 


216 


HD 112999 


96 


102. 


140. 


231 


HD 114886 


92 


109. 


151. 


214 


HD 115071 


71 


95. 


133. 


208 


HD 115455 


81 


96. 


128. 


200 


HD 116852 


66 


98. 


147. 


200 


HD 122879 


90 


105. 


148. 


200 


HD 124314 


74 


98. 


138. 


208 


HD 137595 


72 


94. 


124. 


197 


HD 140037 










HD 144965 


70 


91. 


125. 


203 


HD 147683=^ 


58 


85. 


116. 


185 


HD 147888 


44 


90. 


110. 


181 


HD 148937 


69 


97. 


132. 


228 


HD 152590 


64 


87. 


125. 


205 


HD 152723 


76 


96. 


141. 


201 


HD 154368 


47 


95. 






HD 157857 


86 


99. 


133. 


203 


HD 163758 


79 


142. 


204. 


277 


HD 177989 


49 


85. 


127. 


198 


HD 190918 


102 


156. 


214. 


310 


HD 192035 


68 


92. 


126. 


205 


HD 195965 


91 


103. 


136. 


214 


HD 198781 


65 


92. 


128. 


191 


HD 200775 


44 


104. 


104. 


168 


HD 203532 


47 


78. 


102. 


169 


HD 208905 


77 


97. 


132. 


214 


HD 209481 


78 


97. 


137. 


215 


HD 209975 


73 


104. 


149. 


243 


HD 210121 


51 


83. 


108. 


178 


HD 210809 


87 


126. 


187. 


278 


HD 220057 


65 


87. 


122. 


192 


HD 303308 


91 


121. 


177. 


300 


HD 308813 


73 


92. 


129. 


181 



2.7 
3.3 
2.7 
3.0 
3.1 
3.4 
4.0 
4.4 
3.1 
3.4 
3.1 
3.4 

5.3 
3.8 

2.7 
3.2 
2.7 
4.1 
2.9 
3.6 
2.7 
2.8 
3.3 
3.1 
3.1 
3.8 
3.7 
3.7 
3.6 
2.7 
3.2 

3.0 
2.7 
3.0 
3.1 
3.7 
2.9 
3.2 
2.9 
3.8 
3.9 
2.9 
4.3 
5.2 
13.6 
3.7 
4.1 
4.0 
3.0 
4.6 
4.0 
3.3 
2.7 
3.2 
3.0 
3.4 

5.3 
6.0 
2.9 
2.9 

3.1 
3.0 
3.1 
3.8 



6.9=^ 
9.0 



HD 147683 presents higher- J hnos with To4(CO) = 7.7 and To5(CO) = 8.5 K (see Fig. 1). 
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TABLE 7 
Chemical Results from CN Chemistry 



Star 


Cloud^ 


Afo(CH) 


No{C2) 


A^p(C2) 
(10'^ cm-^, 


iVo(CN) 


Afp(CN) 


T 
(K) 


mv 


nH(CN) 
(cm-») 


BD+48 3437 




4.80 




4.5 


0.58 


0.58 


65 


2.17 


475 


BD+53 2820 


+0.9 


1.59 




<3.9 


<0.93 


<0.93 


65 


2.17 


<1600 




+6.9 


2.49 




<5.6 


<0.93 


<0.93 






<1425 


HD 12323 


-13.5 


1.77 




4.1 


0.89 


0.72 


65 


2.23 


~1600 




-9.7 


3.36 




6.3 


1.03 


1.03 






1200 


HD 13268 


-16.4 


1.20 




<2.1 


<0.42 


<0.42 


65 


2.48 


<1450 




-10.4 


2.70 




3.5 


0.61 


0.60 






900 




-7.4 


5.42 




7.5 


1.33 


1.35 






1000 




-1.0 


1.82 




<2.4 


<0.42 


<0.42 






<925 


HD 13745 


-43.9 


5.89 




6.5 


0.74 


0.74 


65 


1.34 


850 




-18.1 


3.98 




<4.9 


<0.57 


<0.56 






<950 


HD 14434 




9.32 




2.7 


0.55 


0.58 


65 


3.72 


125 


HD 15137 


-13.4 


2.24 




<2.0 


<0.42 


<0.41 


65 


2.98 


<425 




-7.4 


1.62 




<1.8 


<0.42 


<0.42 






<900 




-0.2 


2.56 




<2.1 


<0.42 


<0.42 






<550 


HD 22951/40 Per 




12.0 


3.6 


4.4 


0.64 


0.42 


40 


1.49 


225 


HD 23I8O/0 Per 


+4.6 


7.01 


<4.0 


<3.3 


<0.30 


<0.32 


40 


1.86 


<200 




+7.3 


11.96 


23.0 


16.0 


1.65 


1.91 






625 


HD 23478 


+4.1 


13.53 


7.8 


8.4 


1.03 


0.85 


50 


1.67 


325 




+7.7 


4.72 


6.2 


6.6 


0.80 


0.74 






775 


HD 24190*^ 




9.5 




7.3 


0.75 


0.73 


40 


1.74 


375 


HD 24398/C Per 




22.0 


35.0 


35.0 


3.9 


3.9 


30 


2.05 


700 


HD 24534/X Per 


+5.0 


6.0 


3.4 


3.5 


0.70 


0.62 


20 


3.84 


250 




+7.1 


25.9 


31.0 


31.0 


6.6 


6.7 






650 


HD 24912/5 Per'' 




12.0 


7.9 


4.6 


0.26 


0.49 


70 


1.44 


~250 


HD 27778/62 Tau 


+5.0 


12.8 


24.0 


36.0 


8.9 


5.8 


50 


2.29 


~1100 




+7.2 


9.4 


14.0 


30.0 


5.1 


3.6 






~575 


HD 30122 




15.72 




10.5 


1.58 


1.54 


65 


2.48 


400 






15.72 




6.6 


1.58 


1.53 




3.72 


200 






15.72 




13.5 


1.58 


1.62 




2.11 


400 


HD 36841 




9.92 




7.2 


1.56 


1.54 


65 


3.16 


475 


HD 37367 


+3.8 


3.23 




<2.3 


<0.34 


<0.34 


65 


2.48 


<425 




+6.2 


9.76 




<2.7 


<0.34 


<0.35 






<150 


HD 43818/11 Gem 


-7.0 


3.80 




<2.6 


<0.61 


<0.61 


65 


3.41 


<425 




-3.9 


2.69 




<2.4 


<0.61 


<0.61 






<625 




+1.2 


2.79 




<2.4 


<0.61 


<0.61 






<600 




+5.2 


2.02 




<2.2 


<0.61 


<0.61 






<900 


HD 58510 




5.11 




<4.6 


<0.51 


<0.50 


65 


1.92 


<375 


HD 63005 


+14.3 


4.64 




10.5 


1.40 


1.41 


65 


1.74 


1300 




+21.0 


4.37 




<8.1 


<1.03 


<1.03 






<1025 


HD 96675 




22.76 




49.7 


6.26 


6.25 


50 


1.58 


1425 


HD 99872 




12.6 




<6.1 


<0.62 


<0.63 


65 


1.90 


<200 


HD 102065 


+1.0 


1.1 




<1.4 


<0.62 


<0.17 


65 


0.90 


>1600 




+3.8 


6.0 




<7.9 


<0.94 


<0.93 






<1600 


HD 115455 




17.0 




<8.1 


<1.6 


<1.6 


65 


3.16 


<275 


HD 137595 




12.2 




<7.7 


<0.80 


<0.82 


65 


1.49 


<400 


HD 148184/x Oph 




34.0 


35.0 


19.0 


1.3 


2.7 


60 


2.30 


~300 


HD 149757/C Oph'' 




25.0 


18.0 


21.0 


2.6 


2.2 


60 


1.98 


325 


HD 154368 


-13.1 


2.1 




0.53 


0.21 


0.20 


50 


4.77 


90 




+3.3 


54.1 


51.0 


58.0 


27.0 


22.0 






750 


HD 157857 


+0.0 


5.02 




<6.0 


<1.15 


<1.16 


65 


2.67 


<900 




+4.2 


2.83 




<4.9 


<1.15 


<1.10 






<1600 


HD 185418'' 




13.0 


<10.0 


<1.5 


<0.50 


<0.48 


65 


4.46 


<30 


HD 190918 


+2.1 


1.73 




<1.9 


<0.34 


<0.34 


65 


2.54 


<775 




+18.3 


1.16 




<1.8 


<0.34 


<0.34 






<1175 


HD 192035 


+1.4 


3.31 




<4.3 


<0.53 


<0.53 


65 


2.05 


<575 




+5.6 


11.42 




33.4 


5.07 


5.11 






1550 




+9.4 


2.62 




<4.1 


<0.53 


<0.53 






<725 


HD 192639'' 




28.0 


<10.0 


<3.2 


<0.70 


<0.76 


65 


3.97 


<40 


HD 198781 




13.19 




13.1 


3.29 


3.25 


65 


3.26 


750 


HD 208440'' 




11.7 




<8.2 


<0.90 


<0.87 


65 


1.80 


<325 


HD 208905'' 




5.4 




<7.6 


<0.90 


<0.90 


65 


1.60 


<850 


HD 209399'' 




7.9 




<8.1 


<0.90 


<0.91 


65 


1.28 


<825 


HD 209975/19 Cep'' 




8.5 




<8.1 


<0.90 


<0.90 


65 


1.52 


<600 


HD 210121 




286.2 


65. ( 


) 55.1 


17.35 


19.10 


50 


3.35 


1425 


HD 210809 


-0.6 


1.61 




<4.4 


<0.64 


<0.64 


65 


1.98 


<1325 




+3.4 


3.86 




<5.4 


<0.64 


<0.64 






<575 


HD 217035A'' 




16.8 




<9.0 


<0.90 


<0.90 


65 


2.47 


<125 


HD 220057 


-1.8 


9.33 




11.8 


1.36 


1.37 


65 


1.61 


750 




+1.9 


3.86 




<5.4 


<0.65 


<0.65 






<850 



Note. — All calculations employ Juv ^ 1, except for HD 27778/62 Tau and HD 210121, where /uv ^ 0.5. 
^ If more than one cloud containing CN appear along a line of sight, they are designated by fLSR values from Table 
3, having identical values for T and tuv- Results are for line of sight because some input data are not available 
for all components. 



28 



SheflFer et al. 



TABLE 8 
Chemical Results from CH+ Chemistry 



Star'' 


Af(CH+) 


N{CR) 


/(Ha) 


mv 


nH(CH+) 




(1012 cm-2) 


(1012 cm-2) 






(cm-3) 


BD+53 2820 


6.40 


4.08 


0.11 


2.17 


2.1 


CPD-69 1743 


15.0 


2.20 


0.13 


1.30 


1.0 


CPD-59 2603 


16.0 


12.0 


0.09 


2.42 


2.4 


HD 15137 


14.0 


6.42 


0.24 


2.98 


0.3 


HD 30614/a Cam 


20.0 


6.80 


0.36 


1.98 


0.4 


HD 37367 


32.0 


13.0 


0.30 


2.48 


0.4 


HD 58510 


12.0 


5.11 


0.15 


1.92 


1.4 


HD 93840 


3.30 


1.80 


0.03 


0.52 


31.4 


HD 99872 


23.0 


12.6 


0.28 


1.90 


1.0 


HD 102065 


11.0 


7.10 


0.68 


0.90 


1.3 


HD 114886 


19.0 


9.00 


0.19 


2.48 


0.7 


HD 115455 


17.0 


17.0 


0.23 


3.16 


0.6 


HD 116852 


3.50 


1.70 


0.13 


1.36 


3.1 


HD 122879 


12.0 


3.10 


0.20 


2.17 


0.5 


HD 124314(-19) 


8.60 


1.80 


0.33 


1.37 


0.5 


HD 124314(+2) 


6.20 


7.10 


0.22 


1.37 


4.2 


HD 137595 


18.0 


12.0 


0.60 


1.49 


0.8 


HD 145502/1/ Sco 


6.30 


5.90 


0.10 


1.17 


9.4 


HD 149038/At Nor 


35.0 


10.0 


0.36 


2.36 


0.2 


HD 152590 


19.0 


10.0 


0.22 


2.67 


0.5 


HD 152723 


11.0 


6.40 


0.13 


2.60 


1.1 


HD 157857 


20.0 


7.85 


0.33 


2.67 


0.3 


HD 163758 


1.40 


2.20 


0.08 


3.07 


3.4 


HD 164353/67 Oph 


7.40 


4.50 


0.26 


0.74 


3.5 


HD 185418 


12.0 


13.0 


0.40 


4.46 


0.1 


HD 190918 


14.0 


2.89 


0.07 


2.54 


0.8 


HD 192639 


41.0 


28.0 


0.34 


3.97 


0.1 


HD 203532 


3.00 


25.0 


0.34 


1.92 


11.5 


HD 208440 


8.70 


11.7 


0.21 


1.80 


3.5 


HD 209339 


6.60 


7.90 


0.17 


1.28 


6.5 


HD 209975/19 Cep 


24.0 


8.50 


0.17 


1.52 


1.5 


HD 210809 


7.50 


5.47 


0.10 


1.98 


3.2 


HD 217035A 


21.0 


16.8 


0.38 


2.47 


0.6 


HD 218376/1 Cas 


11.0 


7.60 


0.24 


1.36 


2.4 



^ If more than one eloud containing CN appear along a line of sight, the velocity is 
given in parentheses. 
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0.0 




HD 147683 (15.95; 

I , i_ 



1477.4 1477.5 1477.6 1477.7 1477.8 1477.9 1478.0 

Wavelength (A) 

Fig. 1.— Sample CO spectra (dots) from HST/STIS and Ismod.f fits (soUd lines) are shown for HD 147683 (log N = 15.95) and HD 
122879 (13.11). The second spectrum has been shifted upward by 0.5 continuum units. 



HD 93222 (19.84) 







0.5 



0.0 



1074 



1076 



1078 
Wavelength (A) 



1080 



1082 



Fig. 2.— Sample H2 spectra (dots) from FUSE and Ismod.f fits (soUd lines) are shown for HD 37903 (log N = 20.95) and HD 93222 
(19.84). The second spectrum has been shifted upward by 0.8 continuum units. Each spectrum synthesis with Ismod.f includes also the 
(3—0) and the (4—0) bands. Also seen is one of the CO Rydberg bands at 1076 A. 
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Fig. 3.— Sample optical CH+ spectra (dots) and Ismod.f fits (solid lines) are shown for HD 99872 (from ESO, log N = 13.36) and 
HD 23478 (McDonald, 12.32). The second spectrum has been shifted upward by 0.2 continuum units. Note the presence of two cloud 
components along the sight line toward HD 23478 (also seen in CN in Fig. 5). 
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Fig. 4. — Sample optical CH spectra (dots) and Ismod.f fits (solid lines) are shown for HD 210121 (from McDonald, log N = 13.46) and 
HD 116852 (ESO, 12.23). The second spectrum has been shifted upward by 0.2 continuum units. 



UV Survey of CO and H2 



31 



1.2 



X 

S 1.0 



•' . . * . * 

tr • / • •• 



0) 

> 



0) 



0.8 - 



0.6 



0.4 - 



0.2 



HD 23478 (12.04) 



Rl 



CN B-X (0-0 




HD 210121 (13.20) 



3874.0 



3874.5 
Wavelength (A) 



3875.0 



Fig. 5. — Sample optical CN spectra (dots) from McDonald and Ismod.f fits (solid lines) are shown for HD 210121 (log A'^ = 13.20) and 
HD 23478 (12.04). Both RO and Rl have two cloud components toward HD 23478, as is the case with CH+ in Fig. 3. 
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Fig. 6. — CO versus H2 is shown for our sample of diffuse clouds. Open circles here and in subsequent figures represent sight lines 
probing prominent PDRs. In panel (a) the sample is fit with two power laws, as revealed by t he 10-point means, wi th a break at log N = 
(20.4, 14.1). Panel (b) expands the view to include CO derived for dark clouds (smaller dots. iF^derman et al.|[T99CIIV 
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Fig. 7. — CO versus H2 distrib ution for diffuse and dark clouds is compared with the H (/uv = 0.5), T (1.0), and I (10) theoretical 
models for translucent clouds from Ivan Dishoeck fc Bl ack (1988). Note the overall agreement between the shape of model curves and the 
transition region from diffuse to dark cloud regimes, as well as with the observed slopes of each type of clouds. 



u 

ao 
o 



















1 






• VUV 


detections 










; / 




14 


• Dark 
PDR 


clouds 
targets 












34078 ' •*. j/^ 
, YX 200775 

• • 9^ 


- 
















\*r**' . 


















• 


/• • 


- 












.t 




^m • 
























13 






• 


• 
m 


/ 


-V 


• 


• • 37903 
1 






- 


• . 


^ 


/ 


f 


• 


/ 




- 








--» 1 
















iK* 


m 




• 
• 








- 




• 


A h 






• 








" 




• ; 


4, . 














- 


12 


: ^ . 














±20% 





19 



20 
log N{E,) (cm-^) 



21 



Fig. 8. — CH column density is shown to have a linear correlation with Af(H2), characterized by a single slope of 0.97 ± 0.07. The 
optical-data fit is seen to match the lMattilal II1986I) dark cloud extension of the CH versus H2 relationship. 
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Fig. 9. — CO has a dual-slope relationship with CH as it has with H2 (cf. Fig. 6). The break in slopes is found at log Af(CH) = 13.0 
and log N{CO) = 14.1, in excellent agreement with the CO versus H2 break found in Fig. 6 and with log <CH/H2> = —7.5. 
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Fig. 10. — In panel (a), log Af(CH+) is seen to be correlated with log Af(H2) < 20.3. Similarly in (b), CH+ is seen to be well correlated 
with CO below log N = 14.1. These two power-law breaks are in excellent agreement with the break of CO versus H2 (Fig. 6), as well as 
in excellent mutual agreement that the CH+ abundance stops increasing at <log N> = 13.2 it 0.1. 
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Fig. 11. — Fitting CN detections only versus H2 returns a slope of B = 1.5 it 0.4. (dashed line). A slope of B = 1.8 it 0.3 results from 
using the Buckley-James regression method for censored data. 
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Fig. 12. — Two heavier diatomic molecules, CO and CN, are shown to be well correlated. This confirms earlier clues that CO and CN 
are found together in the same colder and denser clumps of gas. The dashed line shows a fit of CN detections only with B = 1.4 it 0.2, 
and the solid line is a Buckley-James regression that includes all CN upper limits (censored data, B = 1.9 it 0.2). 
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Fig. 13. — Interpolated values of the shielding function of Ivan Dishoeck fc BlacH l|1988l ) are plotted as contours over the observed 
distribution of A'^(CO) versus A'^(H2). UV shielding plays a role in the steepening slope beyond log N{CO) ~ 15. 



^ 



2.6 



2.4 



2.2 



m 

° 2.0 



l.f 







1 








1 ' 


1 




' 1 







a 


To^i^z) 



















_ 


▲ 


^03(^2) 
















_ 


- 





^02(H2) 














D ° 


- 


- 














D 




n 


- 


— 












D 




n 


[IlnD — B— 


-— ' 


- 






n 




a 




-ffn 




iprm a a 


- 


- 




D 








nos^^^ 


U A 


h 


□ 3 


- 


- 


n 




AD 


"^ 


'n 


LP n 

A n 

A A 


u 




A 

AA 
A AA____- 


; 












▲ 


.aA-A 


A 






- 


▲ 


▲ 


""O 




'X' 


■*■ 





— 


_— 




▲ 


▲ 







.^£>-^ — <y 


A 


T' 





- 





^__<>__— 





'^ 


'^ 


^ 0. 


'^0 






- 


- 




1 








1 


1 






1 


- 



1.7 



1.8 1.9 

Log T,,{H,) (K) 



2.0 



Fig. 14. — Three higher- J excitation temperatures of H2 are plotted versus Tex of the J = 1 level, showing three indistinguishable positive 
correlations with < B > = 0.48 ± 0.08. The adopted global uncertainties are ±5% and ±10% in Tqi and Tqj (J > 1), respectively. 
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with either quantity. The relationship has a slope of 0.43 ± 0.06 and r = 0.77. 
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Fig. 16. — Abundance of CO relative to H2 is plotted versus two empirical indicators of gas density: CN/CH in (a) and CN/CH+ in (b). 
The sample of CO/H2 data points shows a tighter correlation with CN/CH+ (r = 0.784) than with CN/CH (r = 0.604). 
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Fig. 17. — Open circles denote values of the density indicator Af(CN)/Af(CH+), showing that higher gas density is associated with 
higher iV(CO) along the upper envelope of the distribution, as well as with higher N(}i2) (from left to right along the diagonal). 
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Fig. 18. — Panel (a) shows density from the CN chemical analysis versus Af(CN). Sight lines with known A''(C2) (filled circles) have r = 
0.61 (solid line) instead of 0.32 (CL increases from 95 to 99%). Panel (b) shows njj versus 7V(CN)/A'^(CH+), where r increases from 0.36 
to 0.64, and CL from 98 to 99.5% owing to restricting the sample to those sight lines with known N(C2)- 
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Fig. 19. — Values of nH(CN) from the chemical analysis are seen in (a) to be well correlated with A'^(C2) values. The slope is B = 0.48 
± 0.15 and r = 0.63. In (b) CN is shown to be linearly related to C2, since the slope of their abundance relationship is 0.97 it 0.28 (r = 
0.64). 




Fig. 20. — Schematic formation routes from C+ to CO involving the most common intermediate gas-phase chemical reactants and 
products. 
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Fig. 21. — CH+ versus H2, CH+ versus CO, and CO versus H2 as a function of nu, for the average value of the far-UV interstellar 
radiation field {I\jy = 1) and AiJ^yri, = 3.3 km s^^. 
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Fig. 22. — CH+ versus H2, CH+ versus CO, and CO versus H2 as a function of /uv> for "H = 100 cm "^ and At^tu^j, = 3.3 km s ^. 
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Fig. 23. — CH+ versus H2, CH+ versus CO, and CO versus H2 as a function of Aiiturb, for /uv = 1 and nj] = 100 cm 



